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I. INTRODUCTION 

RHIC has uncovered an exciting new state of matter, which has partonic degrees of freedom, at ^s„„ ^ 62 to 200 
GeV, the strongly coupled Quark Gluon Plasma (sQGP),. The RHIC-wide consensus on the status of this discovery 
as of a few years ago is documented in the various white-papers from the RHIC experiments [IHl]. In spite of 
impressive progress since 2001, as noted in the white-papers several of the most important questions that motivated 
the construction of RHIC are not yet fully answered, and maybe it is not a coincidence that the available RHIC 
parameter space has not yet been fully explored. 

The QCD phase diagram lies at the heart of what the RHIC Physics Program is all about |5l |B]. While RHIC has 
been operating very successfully at or close to its maximum energy for almost a decade, it has become clear that this 
collider can also be operated at lower energies down to y's^^ — 5 GeV without extensive upgrades. From this purely 
empirical perspective, an exploration of the full region of ^/s~^ available at the RHIC facility is surely imperative. 
The STAR detector, due to its large uniform acceptance and excellent particle identification capabilities, is uniquely 
positioned to carry out this program in depth and detail. The first exploratory beam energy scan (BES) run at RHIC 
took place in 2010 (Run 10), since several STAR upgrades, most importantly a full barrel Time of Flight detector, are 
now completed which add new capabilities important for the interesting physics at BES energies. Our results at top 
energies suggest that a new form of matter, the sQGP, is created and that it is locally equilibrated early-on because 
of its observed hydrodynamic expansion patterns. Also, hadrochemical species equilibrium is observed just after 
hadronization. It appears that the transition to this state, at these high temperatures and low /i^, is a crossover [7], 
i.e. a smooth, continuous transition from a QGP to hadrons. Theoretical model calculations predict that at lower 
temperatures and high baryon chemical potentials this cross-over will become a first order phase transition (8rtl6] 
resulting in a critical point occurring at intermediate temperatures and baryon chemical potentials [5]. 

The data taken via a RHIC BES will be used to explore several of the open questions in the field heavy-ion physics. 
Among the questions are: 

1. Can we see evidence of a Critical Point (CP)? 

2. Can we see evidence of a phase transition? 

3. What is the evolution with ^s„„ of the medium that we produce? i.e. how do the results that indicate the 
presence of the sQGP turn off as ^s„„ is reduced? 

At the forefront of this list is the search for evidence of a CP and/or its associated first order phase transition in the 
phase diagram of nuclear matter. Fig [T] Of course, for any discussion about a phase diagram to be valid, we also 
have to conclusively answer another of the above questions: whether collisions at RHIC form a thermodynamic state? 
While recent progress in lattice QCD and model calculations is indeed impressive, the location of phase boundaries 
between hadronic gas and the sQGP and the exact position, in T and fiB^ of the hypothesized critical point remain 
unknown |17H19| . It therefore falls upon the experiments to find observational evidence of its existence. However, 
the available theoretical estimates indicate that the critical point might be in the region of the phase diagram probed 
by heavy ion experiments (see for example |20H23j ). and in particular within the collision energy range probable 
at RHIC. The proposed Beam Energy Scan (BES) from ^/s^ =5-200 GeV at RHIC is motivated to a considerable 
extent by this exciting possibility of uncovering evidence of a critical point and/or its associated first order phase 
transition line. [248) 

Theory has predicted several signatures of a first order phase transition and for the CP (see for example [MH^U] and 
references therein). Most of these result in increased fluctuations when the freeze-out trajectory passes through the 
CP or strong variations as a function of ^s„„ if the collision energies sampled encompass the range in fiB where the 
CP occurs [Ml EI]. However, the magnitude of these oscillations and the probability of their survival through the 
final re-scatterings in the hadronic state have yet to be rigorously calculated. There is also the question of hard and 
semi-hard processes obscuring the signals of the CP. Some hydrodynamical calculations suggest that the CP acts as 
an attractor [32]. This means that as long as the thermalized medium produced has initial conditions close to the 
CP, as it evolves with time its trajectory through the T, fiB phase space will be focussed towards the CP. Such a 
focussing effect has already been discussed in the context of liquid-gas nuclear transitions [33 . This attraction means 
that the impact due to the inability of the theoretical calculations to pin-point the CP's exact location is minimized, 
at least experimentally, as we only "have to get close" to the correct collision energy to "land" at the CP. A finite 
extension of the critical domain in the /is, T diagram implies that there are unlikely to be any sharp discontinuities. 
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Baryon Cheinical Potential ;Jq (M«V) 



FIG. 1: A schematic of the phase diagram of nuclear matter. The location of the CP is placed within the RHIC BES range. 
Lattice QCD estimates [17H19j indicate that the CP falls within the interval 250 <^b< 450 MeV. The black closed circles are 
current heavy-ion experimental calculations of the chemical freeze-out temperature, Tch, and /xs based on statistical model fits 
to the measured particle ratios. The yellow curves show the estimated trajectories of the possible collision energies at RHIC. 



The size of the correlation lengths, and thus fluctuations, are also restricted, due to the finite system size effects, to 
< 6 fm. A phenomenon known as "critical slowing down" is predicted to cause the correlation lengths to be at most 
2 fm |M1 inn . Despite these problems in precise predictions it is imperative that RHIC lead the way in attempting to 
find the CP. 

Establishing the validity of the CP prediction, or even bounding the region where it is sited by proving the existence 
of both a cross-over AND a first/second order transition, would place RHIC results into text books. Such a result 
would be as seminal as proving the formation of the sQGP. 

In the following sections we describe: 

1. Current proposed measurements, with estimations of the accuracy of the measurements given an assumed event 
count at each yTs^. These measurements include novel extensions to those already made at the SPS. These 
are only possible due to the improved coverage, particle identification (PID) and proposed improved statistics 
over the previous experiments. 

2. The preliminary results from analysis of the data from the RHIC low energy test runs. 

3. The available data to use as the p4-p baselines. 

4. The STAR detector including newly available upgrades pertinent for the energy scan. This section includes de- 
tailed discussion of our particle identification abilities, and our particle acceptance and reconstruction efficiencies 
as a function of multiplicity and pr- 

5. A summary including our proposed running plan 

We conclude with a run plan for data taking that would allow us to accomplish all our initial goals, including 
those mentioned above which are not directly related to the CP search. Should evidence of the CP, first order 
phase transition, or other unexpected results be observed we would then propose further dedicated running at, and 
bracketing, the relevant collision energies. 
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II. PROPOSED PHYSICS MEASUREMENTS 



The physics to be extracted from the BES have been spht into five subsections. First we focus on those measurements 
that seem most hkely to provide evidence of the Critical Point. These are predominantly fluctuation measures. The 
second subsection covers analyses designed to instead provide evidence for a first order phase transition. The signatures 
of the sQGP are covered in sub-section three followed by a sub-section on particle production studies designed, in 
the main, to probe if the medium produced in these collisions is indeed in thermal equilibrium. In the fifth and final 
subsection the study of a potentially sensitive measure of local strong parity violation is presented. 

When discussing performing an energy scan at RHIC the question often arises: why embark on such a scan when 
data has already been taken at the SPS? This section explains how the scan at RHIC will improve on the data taken 
at CERN. Firstly such a program will enable measurements to be made from the SPS energy region on up to the top 
RHIC energies with the same detector. Secondly, since this is a collider experiment the same uniform acceptance will 
occur at each energy point. Finally, not only will the statistics be better due to the higher acceptance of the STAR 
detector, but also, there will be cleaner, and more extensive PID capabilities. This allows us to not only repeat those 
measurements performed at the SPS in much finer detail but also to enhance the studies by performing numerous 
differential measures. It is most likely only by looking differentially that the signals produced by passing close to 
the CP will be extracted. When one is forced to integrate measures, due, for instance to small statistics or limited 
acceptance, much information is lost. 

Many of the potential signals proposed to be resulting from the CP and/or a first order phase transition revolve 
around fluctuation measures. Discussed below are the measures we currently propose. Of course there is still the 
possibility of performing other studies that emerge from more details discussions and continued development of the 
theory. 



A. Locating the Critical Point 



1. Fluctuation Measures 



The characteristic signature of the existence of a CP is an increase of fluctuations |36i- For instance. Lattice QCD 
calculations [37j indicate large fluctuations in the derivatives of the partition functions with respect to baryon, charge, 
and strangeness chemical potentials as a function of the temperature of the system. Of particular interest are the 
moments of the charge/baryon number/strangeness fluctuations which are obtained theoretically from the second 
(quadratic) and fourth (quartic) derivatives of the logarithm of the QCD partition function (x2 and Xi respectively). 
Fig. [2] and Fig. J3] These fluctuations can be related to event-by-event moments of various observables in heavy-ion 
collisions. Fig. [2]and Fig. [3] are calculated for iib ~ 0, such fluctuations are expected to diverge at the Critical Point. 
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FIG. 2: Quadratic fluctuations of baryon number, electric 
charge and strangeness. All quantities have been normal- 
ized to the corresponding free quark gas values and are for 
liB=0 and Tc=200 MeV [37]. 



FIG. 3: Quartic fluctuations of baryon number, electric 
charge and strangeness. All quantities have been normal- 
ized to the corresponding free quark gas values and are for 
Hb^O and Tc=200 MeV 



Measures of particular interest in STAR are fluctuations in the (pt), the K/tt, p/tt, and K/p ratios, and V2. Also of 
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interest are the high moments of the net-protons. 



(Pt) Fluctuations 

(Pt) fluctuations are chaflenging as there are a number of effects that can swamp the signal. For instance eUiptic flow 
can cause a non-statistical fluctuation of the {px) if the experiment does not have 2tt acceptance. Since the plane of 
the collision varies event to event if the acceptance is limited one is forced to measure at a random angle to the event 
plane each event. The {px) in the plane is expected larger than the {px) out of the plane, therefore while the pr 
fluctuations in and out of the plane are independently small the average of a random angle creates an artiflcially large 
apparent fluctuation. This is illustrated in Fig. [4j and shows that for mid-peripheral collisions fluctuations become 
very signiflcant for Au-Au collisions at ^/s~^ =200 GeV. 
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FIG. 4: The (pt) fluctuations as a function of Npart due 
to the presence of elliptic flow. The data are from Au-|-Au 
coUisions at ^/s^ =200 GeV. From [55] . 



FIG. 5: The measured event-by-event (pr) for data (red 
solid circles) , and from event mixing (blue open circles) the 
flt to the data is shown as the red dashed curve, while that 
for the mixed events is the blue solid curve. The data is 
from y/s^ =200 GeV Au+Au collisions. From [39] . 



The (pt) is measured for all events and also estimated for mixed events. The results are compared and any difference 
between the data and mixed events is an indication of non-statistical fluctuations, Fig. [Sj Current results show 
significant non-statistical fluctuations, as measured via the covariance of the two particle transverse momentum 
correlation measure {ApT,i^PT,j) at all energies. Fig. [6]|39|. They increase with .^s„„ and are larger than those 
predicted by HIJING (open symbols in Fig. [6|. The data of Fig. [6] show an initial rapid rise as a function of centrality 
before plateauing around Npart ~ 150, with a common turning point for all measured energies. Again HIJING fails 
to reproduce this trend being essentially centrality independent. If one scales the fluctuations by the (pt) a different 
picture emerges. Fig. [7] The energy dependence is removed and the result is now inversely proportional to the 
centrality. HIJING now reproduces the centrality trend but still underestimates the magnitude. The CERES data 
confirm the energy independence of this measure as shown for the most central data in the inset of Fig. [t] [33] • 

Fig. [Sl shows the RMS of the px fluctuations as a function of .ys„„ . The RMS grows smoothly with increasing 
log(/7sj^j^) [3^. However, there is a relatively large gap in the measurements between 20 and 60 GeV and this region 
ought to be covered before it is claimed that there are no anomalies in this measure. 

An alternative measure of px fluctuations is $(pt) l42j. The difference between event-by-event fluctuations of the data 
and those from mixed events is quantified by ^{px)- If the system created in the collisions does not have inter-particle 
momentum correlations and thus all particles are emitted independently $(pt) will be zero. If instead A-|-A collisions 
are an incoherent superposition of N-|-N events ^{px) will be non-zero but a constant for all centralities in Au-|-Au 
and for p-l-p. For the most central data, once event-by-event impact parameter fluctuations have been removed NA49 
observe no collision energy dependence of ^{px) and its magnitude is consistent with zero 43 . However, they observe 
a signiflcant non-monotonic evolution of $(pt) as a function of centrality for highest SPS collision energies (^s„„ 
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= 17.3GeV as a function of the mean number of wounded 
nucleons (Nw)- Corrections for limited two track resolu- 
tion have been applied. Errors are statistical only. Sys- 
tematic errors are < 1.6 MeV/c. From [41j . 



=17.3 GeV), Fig.|9][lT], this result has been confirmed by CERES [44. Further more detailed studies are needed to 
confirm these results and their potential implications. These will be undertaken by STAR via the BES. 



K/n Fluctuations 



Current STAR results for K/n fluctuations from Au-|-Au collisions at ^s„„ = 19.6, 62.4, 130, and 200 GeV are shown 
in Fig. 10 along with results observed by NA49 at the SPS in central Pb-|-Pb coUisions at ^s„„ = 6.3, 7.6, 8.8, 12.3 
and 17.3 GeV |45]. The fluctuations are analyzed using cTdyr. 
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where adata is 



the relative width (standard deviation divided by the mean) of the K/tt distribution for the data and (Jmixed is the 
relative width of the K/tt distribution for mixed events. Our results for K/tt fluctuations in central collisions show 
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little dependence on the incident energies studied and are on the same order as the NA49 measurements at 
12.3 and 17.3 GeV. 
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FIG. 10: Experimental results for (Jdyn for K/n as a func- 
tion of ^ySj^j^ |45l 146) . Also shown are results from the 
statistical hadronization model of Torrieri [47i . 




FIG. 11: Comparison of the predictions of the HSD and 
UrQMD models to the experimental data for adyn for K/n. 
Data from \M \M 



In Fig. [Toj we compare the statistical hadronization model results of Torrieri [47j to the experimental data. We see that 
when the light quark phase space occupancy, jq, is one, corresponding to equilibrium, the calculations underestimate 
the experimental results at all energies. When 7^ is varied to reproduce the excitation function of /■k'^ yield 
ratios and the excitation function of temperature versus chemical equilibrium over the SPS and RHIC energy ranges 
[171 H5] , the statistical hadronization model correctly predicts the dynamical fluctuations at the higher energies but 
under-predicts the NA49 data at the lower energies, supporting the conclusion that the lower energy fluctuation data 
are anomalous [H]. In these fits, 7, > 1 (chemically over-saturated) for ^/snn < 9 GeV and 7, < 1 (chemically 
under-saturated) for ^/snn > 9 GeV. 

The changes in susceptibilities as a function of temperature, illustrated in Figs. [2] and |3] for /is = 0, are expected 
to diverge at the critical point. Such an effect might be observable as deviations of fluctuations from a monotonic 
dependence on incident energy in central collisions. However, changes in the underlying physics can also induce 
changes in the fluctuations as a function of incident energy. To gain insight into what we might expect from K/tt 
fluctuations as a function of energy, we compare the experimental results to predictions from the HSD model [49] 



and the UrQMD model [SQi in Fig. [llj The NA49 UrQMD resuhs [45j were carried out using UrQMD version 1.3 
with an NA49 acceptance filter while the STAR UrQMD results were carried out using UrQMD version 2.3 with a 
STAR acceptance filter. We can see that UrQMD reproduces the results at RHIC energies, but under-predicts the 
fluctuations at low incident energies. HSD seems to reproduce well the general shape of the measurements, but is 
slightly above the data in the ^s„„ ^ 10 GeV range. The fact that no model completely reproduces the measurements 
at all ysj^„ , combined with the lack of experimental data in the range of ^/s^ =20 - 60 GeV, means that the question 
of non-monotonic behavior of K/n fluctuations must be answered with additional measurements. 

To make these measurements, one needs to attempt to measure all the kaon and pious. The K reconstruction efficiency 
as a function of pr is rather low, see section [VP] and Fig. [77] for more details, due to multiple effects. The first is the 
decay of the K, the dominant channel being — > /Lt^ -I- with a cr of 3.7 m. This means a significant number of 
the kaons decay before being measured by the TPC. Secondly the cuts required to select clean kaons, essential for the 
fluctuation measures, reduce the efficiency further. By using the ToF we can extend the clean PID range to higher pT 
and thus gain essential coverage. As already mentioned not only coverage but clean PID is needed. Again the ToF 
is useful because it can eliminate much of the electron contamination at low pt that dE/dx measures alone cannot 
remove. The most damaging misidentiflcation is that of a X as a tt or vice-versa, since this effect distorts both the 



numerator and denominator of the measure, i.e. K/tt — >■ (K + 1)/{tt— 1) or {K —1)/{tt + 1). Fig. 12 shows a simulation 
of the raw STAR charge integrated K/tt ratio as measured from 100k HIJING simulated central Au+Au collisions 
at y/s~^ =8.8 GeV. The red curve shows the distribution of the K/tt ratio ratio without using the ToF information 
and the black curve shows the distribution of the K/tt ratio using the ToF. It can be seen that the ToF significantly 
improves our ability to study K/tt fluctuations. It has been estimated that 1 M events at each of the proposed collision 
energies for the BES are sufficient to perform these PID fluctuation studies. Figure [T3| shows the estimated statistical 
error for STAR's adyn for the charge integrated K/tt ratio with and without the ToF information. Also shown for 
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comparison are the current NA49 and STAR measured data points. With the ToF, STAR's relative error is ±5%, 
without the ToF, this doubles to ±10%. Although K/tt fluctuations have been studied by NA49 and STAR already, 
there still remains substantial work to fully understand the results. For instance, NA49 only measured the charge 
integrated K/tt fluctuations in central collisions. STAR, with its larger acceptance and ToF identification reach, will 
be able to repeat these measures in much greater detail and also measure the charge separated fluctuations. These 
improvements will allow us to observe event-by-event if there is truly something special happening in the /tt~^ 



ratio as suggested by the integrated measure shown in Fig. 10 



p/tt Fluctuations 

The study of p/tt fluctuations may provide information about baryon fluctuations, p/n fluctuations have been studied 
as a function of y/s~^ by NA49 [33] and by STAR as shown in Fig. fT4| using the variable ddyn at the same energies 
as those used to study K/t: fluctuations. The dependence on ^s„„ oFthe NA49 results and the STAR results seem 
to join smoothly. Here the fluctuations are negative, which may indicate that the decay of resonances is important 



for p/n fluctuations. In Fig. 14 we compare the predictions of UrQMD to the measured fluctuations. The UrQMD 
calculations were done by NA49 45j using the NA49 acceptance filter. We can see that the predictions of the model 
is reasonably close to the experimental results. 

Studying charge separated p/n fluctuations will be crucial to understanding the energy dependence for several reasons. 
One reason is that the p/p changes strongly with ^^/s~/^. Another is that several resonances exist such as A — )• p + n^ 
or A° p + n^ that can have a signiflcant effect on the observed p/n fluctuations. It has been estimated that 1 M 
events at each of the proposed collision energies for the BES are sufficient to perform these PID fluctuation studies. 



K/p fluctuations 

The correlation between strangeness S and baryon number B is sensitive to the state of matter created in heavy-ion 
collisions [53]. In a system consisting of quarks and gluons, strangeness is carried by s and s quarks, which also carry 
fractional baryon number ( B = 1/3 for s and -1/3 for s quarks). In a hadron gas, much of the strangeness is carried 
by kaons for which B = leading to a weaker correlation between strangeness and baryon number, which has been 
reproduced in the string hadronic model UrQMD [Mj. At the transition from hadron gas to quark gluon-plasma, a 
rapid fluctuation in the baryon-strangeness is expected, and recent lattice QCD studies [37| conflrm this prediction. 
As both baryon number and strangeness are conserved quantities, the effect is conveyed to the flnal state and can be 
observed via the fluctuations in the event- wise kaon to proton ratio. The feasibility of measuring K/p fluctuations. 
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(p/K) as a function of 



from NA49 and STAR fSll 



characterized by adym have recently been shown by STAR and NA49, the preliminary results are reproduced in 



Fig. 15 [51] [52]. The same mixing procedure has been used in both of these analyses. There appears to be only a 
slight increase in the dynamical fluctuations at RHIC compared to the top SPS energies. An interesting point to note 
is that adyn for p/K and for K/p are not exactly the same. The low multiplicity of p and K in each event, and that 
o'dyn itself is small, may account for much of the observed differences. More data are needed to shed light on the 
RHIC results. The ToF upgrade will facilitate this important measurement in the beam energy scan. It has been 
estimated that 1 M events at each of the proposed collision energies for the BES are sufficient to perform these PID 
fluctuation studies. 



Net proton Kurtosis 



To date most experimental fluctuation measures have concentrated on the second moments (proportional to the square 
of the correlation length). However, estimates of the magnitude of correlation length in heavy-ion collisions indicate 
that they could be small around the critical point (of the order of 2-3 fm) [S^ , making it challenging to be detected 
in experiments. 

Higher moments of event-by-event pion and proton multiplicities might be significantly more sensitive to the existence 
of the critical point than measures based on second moments. The fourth moment, the kurtosis, of these multiplicity 
distributions is expected to be proportional to the seventh power of the correlation length [SS]. In addition, it is 
expected that the evolution of fluctuations from the critical point to the freeze-out point may lead to a non-Gaussian 
shape in the event-by-event multiplicity distributions. The kurtosis of multiplicity distributions might then provide a 
more sensitive observable for the search of the QCD critical point. 



Further in Lattice calculations, which assume the system to be in thermal equilibrium, the kurtosis of event-by-event 
net-baryon number, net charge and net strangeness are related to the respective susceptibilities. These susceptibilities 
show large values or diverge at the critical temperature [121 [Ml ISZl ES] , Fig. 16 The measurement of higher moments 
of event-by-event identified particle multiplicity distributions and its variation with centrality and beam energy will 
provide the first direct connection between experimental observables and Lattice Gauge Theory calculations |37j . 



Using star's excellent mid-rapidity PID capabilities, (anti) protons for pr =0.2-1.0 GeV/c can be cleanly identified 
using dE/dx in the TPC. We can therefore carry out a proton kurtosis analysis It is advisable to do this analysis 
over as large a px range as possible to preserve the long-range correlations that the kurtosis is most sensitive to probe. 
Therefore the sensitivity of these measurements is greatly increased with the inclusion of the ToF, and measurements 
related to isospin susceptibilities and strangeness susceptibilities will be also possible. 



The UrQMD model [501157] can serve as a baseline of hadronic processes involved in heavy-ion collisions. To estimate 
the minimum statistics required to carry out kurtosis analysis in the RHIC beam energy scan, simulations have been 
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performed in the latest UrQMD version 2.3 [SSUSS]. Figure 17 shows the net-proton kurtosis evaluated within STARs 
central barrel acceptance (|y| < 0.5) for central Au-|-Au collisions in the energy range 5 < ^s„„ < 200 GeV. The 
analysis shows that a measurement of the net proton kurtosis can be made with a statistical error of 0.75 when 100 k 
events are used. Extrapolating these results, we estimate that 4 M events are needed at each beam energy to obtain 
a statistically precise measurement of better than 10%. 



Azimuthal Correlations and V2 Fluctuations 



V2 fluctuations cr^^ have been studied by STAR ^60j and PHOBOS |ST] as a method to test models of the initial 
conditions in heavy-ion collisions and to infer the effectiveness of the conversion of eccentricity to momentum-space 
anisotropy on an event-by-event basis. After further investigation, it has been found that it is not experimentally 
possible to determine a^^ independent of non-flow correlations (S) (see for instance [6 2) and references therein for 
more discussions). The sum of non-flow and elliptic flow correlations (crtot) i'^ the form 5 + '2,a'^^ can be measured via 
the non-statistical width of the flow vector distribution which is related to the difference between the two- and four- 
particle cumulants: a^^^ w S + 2ay^ ~ V2{2}^ — V2{4:}^. Non-flow due to HBT and resonances can be nearly eliminated 
by selecting like charge particles and by using rapidity gaps (particles are required to be well separated in rapidity). 
Those correlations due to jet production should fall off with energy. The energy dependence of a^^^ S + 2(t^^ should 
therefore be of considerable interest in an energy scan as it should be increasingly dominated by a^^ as the beam 
energy is reduced and can be expected to rise if matter is created near the QCD critical point. 



Figure 18 shows preliminary STAR results for af^^ based on 15 million Au+Au collisions at 200 GeV 60 . The 



errors are dominated by systematic uncertainties. Measurements can be carried out using several hundred thousand 
minimum-bias events. This estimate is supported by the analysis of the two- and four-particle cumulants carried out 
on the 130 GeV data based on several hundred thousand events. These measurements will therefore be possible at all 
y^s„„ values of the proposed BES, with likely exception of 5 or 6 GeV. 



Photon Multtphcity Fluctuations 



In a thermodynamical picture of the system formed in the collision, the fluctuations in particle multiplicities can be 
related to the matter compressibility |63l . which could aid our understanding of the critical fluctuations at the QCD 
phase boundary. Through event-by-event photon multiplicity measurements using the Photon Multiplicity Detector 
(PMD) at forward rapidities, we can access, at the same beam energy, multiplicity fluctuations at a higher baryon 
chemical potential compared to those measured at mid-rapidity via charged particles in the TPC. The excellent 
spatial resolution of the PMD allows one to study these fluctuations in localized regions of phase space. Due to the 
uncertainty in the location of the critical point we need to investigate as much of the QCD phase diagram as possible. 
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Given the experimental limitations on the number of beam energy points that can be scanned, the ability to study 
fluctuations across a wide range in rapidity enhances the sensitivity of our search. 



The Fig. 19 shows the relative fluctuation in photon multiplicity {uj^ = a /mean) as a function of number of partici- 
pating nucleons. The results are at a center of mass energy of 17.3 GeV in Pb+Pb collisions within a pseudorapidity 
coverage of 2.9-4.2 measured using the PMD in WA98 experiment |64) . 
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tuations in photon- like clusters and photons after correc- 
tion for efficiency and purity. These are compared to cal- 
culations from a participant model and those from VENUS 
event generator |64) . 



2. Forward-Backward Correlations 



Long-range multiplicity correlations have been measured to provide insight on the mechanisms of particle production. 
Many experiments show strong short-range correlations (SRC) over a region of ± 1 unit in rapidity [65 68 Correlations 
that extend over a longer range are observed in hadron-hadron interactions at higher energies [66, 68j. It has been 
suggested that long-range correlations (LRC) might be enhanced in hadron-nucleus and nucleus-nucleus interactions 
due to multiple parton collisions [SSI [70] . 

The Forward-Backward (FB) correlation strength for multiplicity-multiplicity correlations is defined by the dependence 
of the average charged particle multiplicity in the backward hemisphere, (Ni,), on the event multiplicity in the forward 
hemisphere, {Nf), such that (iVf,) = a + (Nf), where a is a constant and the correlation strength is 



{Nfm)-{Nf){N,) 
(iV|) - {Nf)^ ■ 

For STAR, the forward-backward intervals are located symmetrically about midrapidity {rj = 0). 

Previous studies of FB multiplicity correlations indicate that both short and long range correlations are present 
in Au-|-Au collisions at y^s^^^ = 200 GeV, while only short range correlations are present in the more peripheral 



collisions [7Tk73| . Measurements for 200 GeV Au-|-Au collisions are shown in Fig. 20 along with measurements in 
p-|-p for reference |74| . From the large value of the correlation strength at large Arj, one can infer that a dense partonic 
system is created with multi-parton interactions that lead to long-range correlations in heavy-ion collisions. 
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FIG. 20: FB correlation strength as a function of A77 (a) for Au+Au at the 0-10% most central events, (b) more peripheral 
Au+Au collisions and c) p+p collisions at .^s„„ — 200 GeV. Figure from |74| . 



An important extension of this work will be to examine Au+Au collisions as a function of collision energy, particularly 
in the 5-39 GeV range. By studying the energy dependence of the FB correlation strength, we will 1) be able to test 
models that have been proposed to describe our data (e.g. Parton String Model [75], Color Glass Condensate |76j . 
HIJING ^77j) and 2) search for non-monotonic trends in shape and magnitude of the FB correlations that might signal 
fluctuations due to a critical-point or a first order phase transition at freeze-out. These studies can be carried out 
with approximately 1 million events at each collision energy. 



3. The Focussing Effect of the QCD Cntical Point 
Asakawa et al. |32j proposed that a critical point would act as an attractor of nearby system expansion trajectories 



in the (/xb,T) plane. Fig. 21 illustrates the situation for three system expansion trajectories which originate from 
slightly spaced initial conditions. By considering the different options of a simple crossover transition, a first order 
phase transition, and a similar transition as displaced by an assumed attraction of the trajectory by an adjacent 
critical point, the authors demonstrate that different initial system trajectories may end up in the same hadro-chemical 
freezeout point. 

Different initial system trajectories imply different "starting conditions" , as governed solely by the the collisional center 
of mass energy. One expects that such starting points are determined by the primordial conditions, prevailing during 
the early equilibration phase, just after the end of target-projectile interpenetration. They are thus initially insensitive 
to the particular conditions, encountered later, in the vicinity of the parton-hadron phase transformation line. The 
system expansion trajectories, labelled by their initial entropy to baryon number ratios (as specified in hydrodynamic 
models underlying such considerations), should therefore be shifted smoothly toward lower baryochemical potential, 
with increasing collisional energy. 
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FIG. 21: The QCD phase diagram showing hydrodynamical trajectories with and without a CP. The Solid line - Crossover 
transition, dash-dotted - first-order transition, dashed line with CP. From |32| . 



The proposed study rests on the following assumption: if no special singularities arise along the parton-hadron 
coexistence line (such as a critical point, and/or associated transition from crossover to a first order nature of the 
phase transformation) the smooth spacing with energy of the initial system expansion trajectories should in turn 
lead to a similar smooth spacing of the hadronic chemical freeze-out points. These points in the (/iB,T) plane 
are obtained from the statistical, grand canonical hadronization model [78H82] which analyzes the observed sets of 
hadronic production ratios, at various energies. This analysis derives, for each energy, a " hadrochemical freeze-out 
point" in the (^b,T) plane. These points are shown in Fig. [ij The expanding system must to go through these 
points at hadronization. Thus these points, in their succession with incident energy, mark an entry concerning the 
position of the respective expansion trajectory in this plane. If trajectory re-focusing occurs, owing to the attractor 
mechanism proposed as a consequence of a critical point, the smooth dependence of the freeze-out points on initial 
collisional energy should be re-shuffled. 

What is required for a search concerning this fascinating conjectured mechanism is a systematic precision measurement 
of the hadron production yields (including the hyperons and anti-hyperons) at each incident energy. At the relevant 
low energies of the intended energy scan, the previous data from the SPS low energy runs do not offer adequate 
event statistics. The resulting, apparently smooth succession of freeze-out points, as illustrated in Fig. [T] need to be 
revisited by high statistics STAR runs at these energies. 



B. Locating a Order Transition and/or Changing the Equation of State 



1. Elliptic Flow 



Elliptic flow (V2) measurements have been performed over three orders of magnitude in y/s~^. The most plentiful 
systematics are available for unidentified particles. Figure 22 shows the least detailed measure of anisotropy - pt- 
and particle- type-integrated V2- over a broad range of energies |83j . The strong non-monotonic behavior observed in 
the excitation function reflects an evolution in the driving physics. At the lowest energies, the colliding nuclei orbit 
each other and may fuse to form a rotating compound nucleus which evaporates particles in-plane (v2 >0) due to 
angular momentum considerations. As the beam energy increases, the crossing timescale becomes much shorter than 
the thermalization time, and the passing beam nucleons are too fast to be captured by the target nucleus' mean field; 
at this point, the system is appropriately described in the language of participants and spectators. At beam energies 
~ 200 AMeV, the passing spectators absorb particles emitted by the participant zone in the reaction plane, leading 
to out-of-plane "squeeze-out" [84], and V2 >0. As the energy is increased to top AGS energy (Ebeam ~10 AGeV) 
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and beyond, Lorentz contraction effects and decreasing crossing timescales gradually render the spectators irrelevant. 
It is at these energies, i/s^^ >5 GeV, that the emission anisotropy reflects that from a hot, compressed anisotropic 
zone, and hydrodynamic expansion is assumed to be the physics driving the evolution. 
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FIG. 22: The charged particle V2 as a function of beam energy. 
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FIG. 23: The charged particle V2/e as a function of the charged 
particle mid-rapidity density. From 84 . 



Above AGS energies any energy dependence is expected to be driven not by changing physics mechanisms, but by 
characteristics - equation of state, viscosity, the number of degrees of freedom etc- of the system determining the 



hydrodynamic evolution. Such effects will be more subtle than the obvious structures in Fig. 22 described above. 
Hence, we will need more detail than integrated V2 values. First, if instead of plotting against ^s„„ one calculates 
V2/e (where e is the initial- state eccentricity) and plots this quantity versus the measured charged particle density 



per unit overlap area, S, one sees that all the data from different energies collapse onto a common curve. Fig. 23 8^. 



Also one can see that it is only at the largest particle densities, corresponding to central collisions at maximum RHIC 
energies, that the data reach the hydrodynamical limit for elliptic flow. This immediately provokes the question about 
the validity of using hydrodynamical calculations for lower collision energies. In light of this and other results much 
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work has been done recently to implement viscous hydrodynamical theories. 

The bulk nature of the hypothesized sQGP phase is revealed in strong elliptic flow, which in central collisions ap- 
proaches the predictions of ideal hydrodynamics, assuming system thermalization on an extremely short timescale 
0.5 fm/c) |85l 186] . However, the mechanism behind such rapid thermalization remains far from clear and is 
under active theoretical study [87l |88j . This has a connection to another novel phenomenon that could be relevant 
at RHIC — saturation of the gluon distribution, often referred to as the Color Glass Condensate or CGC — which 
characterizes the nuclear parton distribution prior to collision ^89 . Various theoretical approaches to connect collision 
geometry, saturated gluon distributions, and the onset of bulk collective behavior are being explored [13 [53110]; more 
experimental input is needed to guide these efforts. 
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If one performs more differential analyses and looks at V2{pt) as a function of ^/s~^ the results shown in Fig 
emerge. For a fixed centrality the V2 at fixed pr values first grows with ^/s^ but then reaches a plateau. This leveling^ 
off at similar y/s^ for differing transverse momenta might present evidence of a softening of the equation of state 
due to a phase transition. It will be extremely enlightening to see if this feature is preserved for differing centralities 
and using identified particle results. Moreover, the energy where the leveling-off begins lies in the unexplored region 
of the RHIC BES, and ought to be determined with more precision. Table [ij shows the number of events estimated to 
be needed in order to measure the inclusive particle V2 up to px ^1.5 GeV/c. 



TABLE I: Estimate of events needed to measure the inclusive particle V2 up to pr ~1.5 GeV/c. 
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FIG. 24: The charged particle V2 as a function of 
for two different pr values. From [83j . 



FIG. 25: Proton V2 as a function of rapidity at ^/s]^ = 
8.77 GeV from Pb-Pb collisions for centrality bin b=3.4-5.3 
fm. From [91]. 



A potential direct signature of a phase transition is the "collapse" of the proton V2 [251 E21 193] • NA49 first reported 
observation of such a result for y/s]^ = 8.77 GeV. However, when they repeated this measurement using other 
techniques to extract V2 the results differ. Fig. 25 [ST]. The results presented in Fig. 25 are from the "standard" V2 
method (blue circles), cumulants for two-particle correlations, V2{2}, (green squares) and cumulants for four-particle 
correlations, V2{4}, (red triangles). The differences between these results are beyond the statistical error bars, however, 
the techniques used are known to have different sensitivities to non-flow effects, such as decays (see [HI] and references 
within for more details). It is currently unclear if these differences are due to non-flow effects or physical fluctuations 
of the V2 that come in with differing weights to these different techniques, or represent an indication of the systematics 
of the measure. With the BES we can perform these measurements with higher statistics and better systematics. Also, 
recent advances in experimental techniques have developed methods to determine the scale of the V2 fluctuations. It 
is therefore essential that we recalculate the V2 and perform the fluctuation measurements with the same apparatus to 
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try and disentangle all these effects. The feasibility of such a measurement is shown in Section III where statistically 
relevant identified V2 measurements have been made with only ~3 k events. 

The first evidence of photon azimuthal anisotropy at SPS energies was observed in the distribution of photons from 
S+Au collisions at SPS energies measured in the pre-shower photon multiplicity detector of the WA93 experiment 
at CERN [94] . Subsequent measurements of both directed and elliptic flow of photons were performed in the WA98 



experiment at SPS, these are shown in Fig. 26 [55]. The preliminary measurements at RHIC using the PMD were 



reported in Ref. |96| . These photon measurements will complement the charged particle measurements at mid-rapidity. 



0.06 r 
0.05 r 
0.04 r 
0.03 r 
0.02 T 
0.01 r 



-A- 



▼ PMD data 
A Monte Carlo 

o V rjt)=o 



^ -t- 



t -t- 



50 100 150 200 250 300 350 400 450 

EJGeV) 



FIG. 26: Directed and elliptic photon anisotropy coeffi- 
cients in the pseudorapidity region 3.25-3.75 for different 
centralities as measured in PMD for WA98 experiment at 
SPS energies [95]. For more details see [95 . 
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FIG. 27: Directed flow of protons from ideal hydrody- 
namics with a QGP phase (open symbols) and from the 
Quark Gluon String Model without QGP phase (full sym- 
bols) [97]. Figure from [98]. 



2. Directed Flow 



Directed flow is obtained from the first harmonic (vi) in the Fourier expansion of the azimuthal anisotropy of emitted 
particles with respect to the collision reaction plane, ^f^, 
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It describes collective sideward motion of the produced particles and nuclear fragments and, even more so than elliptic 
flow, carries information on the very earliest stage of the collision. Specifically, at large rj (in the fragmentation region) , 
directed flow is believed to be generated during the nuclear passage time (2i?/7 ~ 0.1 fm/c) |100| 1101] . It therefore 
probes the onset of bulk collective dynamics during thermalization, providing valuable experimental guidance to 
models of the pre-equilibrium stage. 

The shape of vi vs. rapidity is of special interest because it has been identified in several theoretical papers as a 
promising Quark-Gluon Plasma signature [Ml ri02| . At low relativistic energies, Vi{y) is almost directly propor- 
tional to rapidity, and relative to protons, the pion vi{y) is significantly smaller in magnitude and opposite in sign. 
The sign of vi is conventionally defined as positive by the direction of nucleon flow in the projectile fragmentation 
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region. Often, just the slope of vi{y) at mid-rapidity has been used to define the strength of directed flow. At RHIC 
energies, directed flow is predicted to be small near mid-rapidity with very weak dependence on pseudorapidity. Cal- 
culations involving a QGP phase suggest that vi (y) may exhibit a characteristic "wiggle" [Wl [551 ll02H104j , whereby 
directed flow changes sign three times, not counting a possible sign change near beam rapidities (in contrast to the 



observed sideward deflection pattern at lower energy, where the sign changes only at mid-rapidity), see Fig. 27 if 
a QGP phase transition is assumed. In these calculations, the wiggle structure is interpreted as a consequence of 
the expansion of the system, which is initially tilted with respect to the beam direction. The expansion leads to the 
so-called anti-flow or third flow component. Such a flow can reverse, over a region on either side of mid-rapidity, the 
normal pattern of sideward deflection as seen at lower energy, and hence can result in either a flatness of wi, or a 
wiggle structure if the expansion is strong enough. However, a similar wiggle structure in proton vi is also predicted 
if one assumes strong but incomplete baryon stopping together with strong space-momentum correlations caused by 
transverse radial expansion |103[|104] . The situation for pion directed flow is less clear in such models. While RQMD 
model calculations indicate that shadowing by protons causes the pions to flow mostly with opposite sign to the 
protons, mirroring the proton wiggle, other calculations predict that pions flow with opposite sign only in a limited 
rapidity range. 

It is not until we move down into the collision energy domain of the proposed BES that the acceptance in pseudora- 
pidity of star's Forward TPCs comes close enough to the rapidity region of the incoming nuclei to fully map-out 
all the changes in sign of vi, and therefore, a more thorough understanding of the wiggle phenomenon will then be 
possible. Figure 
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shows the measured directed flow for Au-|-Au collisions at ^s„„ — 9.2, 62.4 and 200 GeV for 
unidentified charged hadrons, and for identified pions in 0-60% centrality Pb-|-Pb collisions at y/s~^ —8.8 GeV |91j . 
For the 62.4 and 200 GeV the data is from 30-60% central data, due to a lack of statistics the 9.2 GeV covers 0-60%. 
This vi study at ^/J/vw = 9.2 GeV, based on just a few hours of sporadic collisions during a machine development 
test in March 2008 has already corroborated the promising nature of this line of investigation. The predicted sign 
change near beam rapidities is clearly evident for the lower energy data but there are insufficient statistics currently 
to attempt to identify any possible "wiggles" at mid-rapidity. 

Based on experimental analyses to date, and on simulations of the performance of the STAR TOF system, a relatively 
modest event sample (0.5 M) — many times smaller than the few million minimum-bias events needed when using 
elliptic flow of abundant baryons and mesons to study constituent quark scaling — will be sufficient to unambiguously 
resolve the pending questions about the "wiggle" phenomenon, discussed above. 



Important insights into the evolution of the system also can be obtained from flow of identifled particles. At low ^s„„ 
most of the protons are transported to mid-rapidity by baryon stopping while pions are created from the collisions, 
and the difference of their flow pattern will shed light on the evolution of the system. In particular, if the vi of protons 
and pions would have the same sign near mid-rapidity, then that would be a signature for a tilted source and would 
be consistent with one-fluid hydrodynamic models incorporating a flrst-order phase transition. It has been argued by 
Stocker that the so-called "collapse of proton flow" reported at CERN may be evidence for such a transition [5H] . All 
of the above considerations stimulate much interest in directed flow over the proposed range of the BES. 

A separate but equally curious sign-change phenomenon occurs for charged-particle Vi{pt) in the vicinity oi px ^ 1 
GeV/c near mid-rapidity at higher RHIC energies (i.e. it is observed in the main TPC of STAR but not in the 
FTPCs). If we assume that vi for pions and protons has opposite sign, then the steeply increasing relative abundance 
of protons above 1 GeV/c has the potential to explain the change of sign, but more data are needed over a wider 
range of beam energies to determine if this phenomenon is properly understood. 



S. Identified Particle pr Spectra 

Changes in shape of the particle pt spectra as a function of ,Js~^ and centrality can be studied to reveal changes in 
the kinetic freeze-out temperature of the hadrons from the medium and their mean transverse velocity. Typically a 
Blast- Wave function is used to extract these variables from flts to the pT spectra |106j . The fits are usually restricted 
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to pt < 2 GeV/c so that the bulk of the particle kinematics are not dominated by hard processes. Fig. 29 shows the 
extracted kinetic freeze-out temperature, T^m, and mean radial velocity, /3f , for identified pions as a function of 
for the most central events 107, 108^. It can be seen that at very low ^/s~^ there is a steep rise in the extracted T^ 
and /3t. For Tkin a plateau appears at higher ^s„„ , although the radial velocity continues to increase. 
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FIG. 28: The charged particle vi as a function of rj for Au+Au collisions at ySj^„ 
identified pions in y^s^^ =8.8 GeV Pb+Pb collisions [91) . 



=9.2, 62.4 and 200 GeV [105], and for 



If instead of tt one looks at , Fig. 30 a more interesting picture emerges. Instead of the extracted Tkin rising 
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is seen for A 



steadily to the topmost collision energies, the measured inverse slope of the pt distribution appears to form a plateau 
for y/s^ ^ 8-12 GeV before rising to significantly higher values for ^s„„ =130 and 200 GeV |109j . Such behavior 
A collisions but not for p+p. If at intermediate ^s„„ the system instead forms a mixed phase region 



the early stage pressure and temperature are predicted to become independent of the energy density [110,- This 
effect creates a step like dependence of the pressure and temperature on the collision energy. This leads, in turn, to 
a weakening of the increase of the inverse slope parameter with ^s„„ as seen in Fig. 
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Measurements need to be 

made to see how far this possible plateau extends as this may signal the transition to a single phase QGP sys tem . 
The K+/7r+ ratio also exhibits interesting behavior around these collision energies (see subsection II B 4 and Fig. 31 1. 
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FIG. 29: The calculated thermal/kinetic freeze-out tem- 
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FIG. 30: The extracted inverse slope parameter of the pt 
distributions of as a function of ^/s~^ for A+A and 
p+p collisions. From |109| . 
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FIG. 31: The K+ /tt+ and K' /tt' ratio as a function of 
■JSj^j^ [111] . Low energy data added to plot from 



FIG. 32: The A/p ratio as a function of ^^s^. From [78] . 
The solid curve is a prediction from the statistical hadron 
gas model described in [112) . 



4- The K/n and A/p Ratios 

One of the most challenging set of results to model from the SPS are the K+/7r+ ratios as a function of ^s„„ , shown 



in Fig, 31 with data from |1091 1113[lll7] . Given this result, and that of the fitted inverse slopes of the K+ (Fig. 30 1, 
it is essential to go back and study this region in greater detail, as it suggests novel physics occurring. Smaller error 
bars will also prove if this peak is merely a statistical fluctuation. With our increased precision we will also be able 
to study this ratio differentially and see if this peak persists as a function of px and centrality. 

Another ratio that is equally hard to explain with statistical hadronization models is the A/p ratio. This ratio 



becomes extremely high at smaller ^s„„ , Fig. 32 [78]. These measurements, especially those at the AGS, were an 
experimental tour-de-force, but lack statistical precision. STAR's large acceptance and high reconstruction efficiency 
of the A will drastically reduce the systematic and statistical errors. It is estimated that 1 Million events at ^s„„ = 
6 GeV will result in several thousand reconstructed A. 

It is possible that the dramatic rise at low ^s„„ is due to p-p annihilation in the baryon-rich medium produced at 
these energies [1181 1119| . However, it remains to be seen if this can account for all of the observed behavior since 
K+/7r+ and A/p both represent s/d and should therefore carry some of the same physics. It is also evident that there 
is a gap in these measurements between ^s„„ =20-60, under our proposed run plan this region will be filled in. 



5. Interferometry 

By studying interferometry as a function of beam energy we can infer the energy density of the medium produced 
at the last re-scattering of the hadrons. The source dimensions, or homogeneity regions, as determined via HBT 
encode different information. Kside only contains information about the spatial extension while Rout holds spatial 
and temporal data. The ratio Rout/Rside can therefore reveal the emission duration of the source. It is predicted 
that for a first order phase transition this ratio should get much greater than unity due to a stalling in the emission 
during the phase transition [120J. Measurements of Ront: R-side and R/ong as a function of collision energy have been 



studied in great detail and are shown, including new results at y/s^ = 9.2 GeV, later in Fig. 61 No major jumps 
in Rotit/Rside are observed. 
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FIG. 33: HBT radii relative to 2nd-order event plane, in Au+Au at 200 GeV [HT] . 
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FIG. 34: Freeze-out anisotropy from 2" -order oscillations of HBT radii. Inset shows 
hydro, evolution of the source shape for an equation of state with (upper) and without 
(lower) softening due to finite latent heat [122j . 
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Momentum spectra and anisotropy tell only half of the story of collective flow. The bulk response of the system 
has a non-trivial structure in both space and time. Just as the pr-dependence of azimuthally-integrated HBT radii 
gives access to the geometric substructure generated by radial flow (e.g. |123j ). HBT measured relative to the first- 
and second-order event plane are the spatial analogs of directed and elliptic flow |1231 1124j , respectively, and contain 
important information not accessible in momentum space alone. As explained below, these measurements can be 
sensitive to a softening in the equation of state, related to a first-order phase transition, or even rapid crossover. 



STAR has measured oscillations of pion HBT radii relative to the second-order event plane (Fig. 33 1 |121j . In addition 
to the overall size of the source, these reveal that the transverse shape is extended out of the reaction plane at freeze- 
out (the stage probed by HBT). It is, however, less anisotropic (more round) than the initial source defined by the 
overlap of the colliding nuclei at finite impact parameter, reflecting the evolution over time of preferential in-plane 
expansion. The anisotropy has been measured at a few lower energies, as well. Since the lifetime of the system and 
the elliptic flow increase with collision energy, one naively expects that, for a fixed initial anisotropy, the freeze-out 
anisotropy becomes less and less out-of-planc extended, and may even become in-plane extended (as predicted |125j 



for example at the LHC). Figure 34 shows the freeze-out anisotropy calculated from the second-order oscillations 
of these HBT radii as a function of collision energy. The results show an intriguing non-monotonic behavior. A 
possible explanation may be as follows: at low energies (say 3-10 GeV) the stiff equation of state of a hadronic system 
generates a large pressure, pushing the system quickly towards a round shape. But at some energy (say 20 GeV), a 
threshold to generate a phase transition is crossed, characterized by a finite latent heat. This generates a soft point 
in the equation of state, and the push towards a round, e = 0, state stalls briefly. As the energy increases beyond this 
threshold, the time spent in the soft state grows, and the out-of-plane-ness at freeze-out grows with energy, until some 
point (say 70 GeV). Then, at even higher energy, the system spends most of its time in the (stiff) QGP phase, and 
the out-of-plane-ness again decreases with energy with no further non-monotonic behaviour. This would be the direct 
analog of the non-monotonic excitation function of V2 originally predicted by ideal hydro models with a softening 
due to a phase transition [93] . The signal in V2 has not been observed, perhaps because increasing viscous effects at 
lower energies smears the structure. However, the spatial anisotropy probed by HBT is weighted in the time evolution 
differently, so may retain sensitivity to the softest point. Figure |34] represents one of the very rare bulk-sector probes 
with a non-monotonic excitation function. Especially given its potential to probe the long-sought "soft point, this 
excitation function must be mapped. 



Measuring HBT correlations relative to the first-order event plane yields even more unique geometrical information. 
In particular, if one approximates the spatial configuration of the freeze-out system as an ellipse, one can extract the 
angle between its major axis and the beam direction. This tilt angle |124j is the spatial analog of the so-called flow 
angle |126| formerly used to characterize directed flow. It is, however, much larger (~40'^ at AGS energies, compared 
for flow angles of < 2'^) and may even have opposite sign. Simultaneous measurement of both the tilt and flow 
angles provides unique insight on the nature and physics behind directed flow at lower energies. At RHIC energies, 
the directed flow becomes even more important, since here we are studying bulk response of the system at the very 
earliest stages of the collision. As discussed above, crossing a threshold to a phase transition will generate a wiggle in 
the directed flow at midrapidity, as the system spends most of the relevant time (keeping in mind that the very small 
time window to generate directed flow is of order of the crossing time) in the soft state. This same physical scenario 
is predicted to generate a non-trivial fingerprint on the coordinate-space configuration (of which the tilt angle is the 
dominant component). The geometry will probe the physics behind the third component of flow generating the vi 
wiggle. Table [H] shows the estimated number of events required to perform these measurements. 



TABLE II: Estimate of number of events needed to measure azimuthally sensitive HBT. 
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11.5 


17.3 
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39 
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C. Turning off the Quark Gluon Plasma Signatures 



RHIC has reported several measures that indicate that the medium created passes through a partonic phase at top 
energies [THl] . We are interested in studying how these signals evolve as a function of y^s^. As well as these partonic 
signatures, there are interesting results as a function of centrality that suggest different physical mechanisms dominate 
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quarks in the hadron as a function of transverse kinetic energy scaled 
by the number of constituent quarks in the hadron for Au+Au colli- 
sions at ,/s7777 =200 GeV. 



when transverse particle densities change from low to high, and the switch from one to the other is both dramatic 
in appearance and sudden. Many studies of these results can be performed with the energy scan. Some of the most 
promising are discussed below. 

1. Constituent Quark Scaling of Elliptic Flow and Elliptic Flow of Multi-Strange Hadrons 

As already stated above, extensive identified particle V2 measurements have been made. Some of these are shown in 
Fig. ?? for various particle species as a function of transverse kinetic energy for Au+Au collisions at y/s~^ = 200 
GeV at RHIC. Baryons and mesons fall onto two separate curves. At low my-mo the V2 follows hydrodynamical type 
behavior |127) . However, at intermediate values of m^-mo the magnitude of the V2 reaches a plateau with baryon V2 > 
meson V2, Fig. ??. If however one scales the V2 and transverse kinetic energy, or px, by the number of constituent 
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quarks in the hadron, all particles now fall on a common curve |128H131] , even in the region where two curves are 
seen in Fig. 36 Such scaling can be explained by quark coalescence, or recombination models |132H135] . which 
provide an intriguing framework for hadronization of bulk partonic matter at RHIC. The essential degrees of freedom 
at the hadronization seem to be effective constituent quarks which have developed a collective elliptic flow during the 
partonic evolution. Further evidence is the sizable magnitude of V2 of the (j) and Q. Both of these hadrons have small 
hadronic cross-sections so it is unlikely that they can develop an elliptic flow during the hadronic phase. Since V2 is 
self quenching this suggests that the early stages of the collision are partonic |1301 1131] . 



If this is indeed a QGP signature, one would expect these effects to turn off at lower y/s~^, especially if we drop below 
the transition temperature at the lowest y/s^. At these energies, the elliptic flow would develop when the degrees 
of freedom of the system are dominated by hadronic interactions. To examine if the number of constituent quark 
scaling would still be observed when in the hadronic stage, we looked at simulations from the AMPT model |136j . 
This model can be run in two modes, with and without a partonic stage. Mode one is without string melting and the 
simulation involves purely hadronic interactions; this is the default setting. In mode two, string melting is included, 
resulting in a fully partonic stage at early times. Full details of the model can be found in reference |136| . AMPT 
with string melting has been shown to reproduce the trends of the RHIC V2 data at y/s~^ — 200 GeV, but it cannot 
when string melting is turned off. The results for identified particle V2 for both these modes at y/s~^ = 9.2 GeV are 
shown in Fig. |37[ Constituent quark scaling is suggested in the string melting scenario whereas it is clearly absent in 
the case where only fragmentation via the Lund string model occurs. 
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FIG. 37: The identified particle V2 parameter from AMPT model at = 9.2 GeV, with and without string melting. See 

text for more details. 



This measurement requires high statistics for identified hadrons at intermediate pT- As the collision energy drops, 
such measurements therefore become very challenging. There are insufficient statistics to make this measurement with 
the current SPS data. Figure [38] shows an estimate of the statistical errors of the V2 for identified protons and tt at 
y/Sj^j^ = 7.7 GeV from 1.4 M 0-80% centrality events from AMPT predictions. The green curves show the estimated 
V2 for baryons and mesons should Number of Constituent Quark (NCQ) scaling exist at ^s„„ =7.7 GeV. Also shown 
in this figure are estimates of the relative error on the proton and tt V2 at pT = 2, 2.5, and 3 GeV/c for 4 M events. 
Extrapolating from this simulations we conclude that for ^/s^^ > 12 GeV and higher a significant measurement of 
potential constituent quark scaling can be made for tt, K (K°), p and, A up to (m7'-mo)/NCQ « 1.5 GeV with ~^5 M 
minimum bias events. For the (j) and ft at least 25 M events are required to reach p-r =2 GeV with a 10% statistical 
error due to their low production rates and the poor signal to background ratios |1371 1138] . 
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FIG. 38: AMPT predictions of identified particle V2 as a function of pr for 1.48 M Au+Au collisions at -^/s]^ =7.7 GeV. 
Also shown is an estimation of the statistical errors for identified proton and tt V2 for 4 M events at 2, 2.5 and 3 GeV/c. The 
solid/dashed green curves are the predictions from NCQ scaling of V2 for baryons/mesons. 



Even if a QGP is always formed at collision energies available at RHIC, the fraction of time spent in the hadronic 
phase should increase with decreasing ^s„„ . If constituent quark scaling is a pure partonic signal this may become 
washed out via hadronic interactions in this later stage. Thus a lack of constituent quark scaling could indicate either 
failure to push the system above the critical temperature or domination by the hadronic stage of the collision on the 
measured signals. 



The PHENIX collaboration has reported an interesting systematic trend in the deviations of the V2 of protons from 
constituent quark scaling as a function of centrality |I39J. STAR have reported similar observations for protons and 
A |140j . To measure the accuracy of the constituent quark scaling with scaled transverse kinetic energy, PHENIX 
perform a fit to all the scaled measured data and then plot the deviations of each hadronic species to this common 
fit. This fit is performed as a function of centrality for ^s„„ = 200 GeV Au+Au collisions and shown in Fig. 39 If 
one focuses on the (anti)proton data (the blue triangles) one sees that for the more central collisions, the low pToata 
fall below the common fit. As the centrality decreases this deviation from the scaling diminishes smoothly until the 
40-50% data is reached at which point the data changes dramatically, to show a positive deviation. The 40-50% data 
is at Npait ~ 90 matching the centrality where the low pr ridge correlation measured by STAR also starts to show 



dramatic deviations from a p+p linear superposition model |141j , and this result is discussed in the sub-section II C 3 
This feature of the low pT identified particle V2 measurements can be extremely well mapped during the energy scan, 
to see if it is a coincidence or whether possible dramatic changes in the underlying physics mechanism affect both 
measures. 



2. Nuclear Modification Factors and Baryon/Meson Ratios 



At intermediate to high pT, it was initially expected that hard processes, that could be calculated via pQCD, would 
dominate. The nuclear modification factor, Rcp, is the '^un scaled ratio of central to peripheral pr distributions. If 
hard processes were not affected by the presence of the medium, and should scale with the number of binary collisions 
and Rcp would equal unity. As shown in Fig. 40 highpT- particles are strongly suppressed |142j . Above 4-5 GeV/c all 
particle species, including non-photonic electrons from the decay of heavy-flavored hadrons, are suppressed by a factor 
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FIG. 39: Identified particle V2 per constituent quark divided by a common fit function as a function of transverse kinetic energy 
per constituent quarks for Au+Au collisions at .ys„„ = 200 GeV for different centralities. 



of 5 in central Au-Au collisions. Attempts have been made to measure the Rcp of hadrons at the SPS |143l 1144] , 
Unfortunately the reach of their data is not sufficient to make a firm statement as to whether the suppression is the 
same as that at top RHIC energies. Since all species measured to date show significant suppression one can measure 
the charged hadron nuclear suppression factor to determine that "jet suppression" has started. At lower collision 
energies initial state effects such as the Cronin effect |145j - the enhancement of particle yield at intermediate pr with 
respect to binary collision scaling - become more prominent. By measuring the nuclear modification factor, Rcp, 
instead of Kaa, the Nbi„ scaled ratio of A+A to p+p pt distributions, the initial state effects should be minimized. 
Finally, since the suppression is also large, one can have a sizable statistical error and still make a meaningful estimate. 
Tabic |III| shows estimates of the number of events needed to determine Rcp of charged hadrons at each beam energy, 
in order to locate the beam energy at which interactions with the medium begin to affect hard partons. It has been 
estimated that too many events are needed to perform this measurement at lower energies. 



TABLE III: Estimates of events needed to measure Rcf up to various pr values at three different beam energies. It has been 
estimated that too many events are needed to perform this measurement at lower energies so they are not shown here. 



Vs^i^ (GeV) 


18 


27 


39 


Pt reach (GeV/c) 


4.5 


5.5 


6.0 


Number of Events 


15 M 


33 M 


24 M 



One of the first indications that particle production at pt above 2 GeV/c was not purely from modified jet fragmenta- 
tion came from the measurement of baryon/meson ratios. The p/tt, and A./K^ ratios surpassed unity for intermediate 
Pt in the more central events, far higher than that observed in p+p collisions, Fig. |4l] [1421 11461 1147j . At higher 
Pt the ratios turn over and decrease steadily until reaching values close to that observed in elementary collisions at 
Pt around 6 GeV/c. This intermediate pt range is also the region where the V2 of baryons is larger than that of 
mesons. Novel hadronization mechanisms are required to reproduce these data. The majority of these models use 
a recombination/coalescence (ReCo) mechanism to form hadrons, as with the quark scaling of V2, described above, 
these models only combine constituent quarks [132H135] . This mechanism naturally leads to an enhanced baryon- 
to-meson ratio when the parton pt distribution is an exponential, i.e. when hadron production is not dominated 
by fragmentation which gives a power-law-like distribution. Such models require quarks to coalesce and hence the 
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FIG. 40: Identified particle Rcp for Au+Au collisions at 



FIG. 41: The A/i^l? ratio as a function of pr for Au+Au 
and p-p collisions at 200 GeV. From |142j . 



observation of ReCo-like behavior is one of the corner-stone pieces of evidence of the formation of the sQGP. 

One of the interesting topics of the energy scan is to investigate in which energies these phenomena are prevalent. 
Although the cross-section for high px processes is smaller at lower y/s]^, the push to high px due to radial flow 
is also lower so the ReCo regime may move to lower pT ranges, counterbalancing somewhat the loss in the hard 
scattering production rate. 



3. Jet Correlations and The Ridge 

Hard scattered partons fragment into a spray of coUimated hadrons known as a "jet" . By examining di-hadron At; 
-A0 (Ary —rii-r]2 and A(t)—4>i-(j)2) correlations, these jets can be identified in the presence of a large background. This 
background is subtracted statistically. In p-p collisions, clear back-to-back peaks from the high interactions are 
apparent. However, in more central A+A collisions, the away-side correlation disappears. This is believed to be due 
to re-scattering of the away-side parton as it traverses the sQGP. The near-side correlation remains, showing evidence 
of vacuum-like fragmentation. While making these measurements, an interesting feature emerged in more central 
A+A collisions, namely a long range A77 correlation at small A(j). This long range correlation sits under the jet peak 
and extends to at least A77=l. This phenomenon, called the ridge, appears to be correlated with jet-like triggers but 
the particles within the correlation have features reminiscent of the bulk rather than jet fragmentation, i.e. the p/tt 
ratio is close to that of the bulk, not a jet, and the Pt spectra are softer than those from fragmentation [HH] . This 
ridge correlation appears in both high p^ triggered and un-triggered correlation studies. 



UnTriggered Correlations 

At RHIC, differential analyses have been developed to allow more detailed investigations of the observed fluctuations. 
Instead of looking at the event-wise quantities, we instead study correlations: 

Ap 



where prcf is the number of mixed event pairs in the distribution being studied and Ap is the number of pairs in 
the real event minus prcf . These correlations are then plotted as a function of relative angles and 77. The resulting 



distributions for Au+Au collisions at -^s^^ = 62 and 200 GeV are shown in Fig. 42 |149j . Clear structures, including 
the near-side jet peak can be seen for all centralities. In p-p collisions, these distributions can be well described by a 
fit function derived from four simple components. These are a: 



1. 1-D Gaussian in 0a accounting for longitudinal fragmentation. 
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FIG. 42: Pair density correlations for Au+Au collisions at ys„„ — 62 and 200 GeV as a function of centrality. From I149| . 



2. 2-D exponential at very small Arj, A(/) for HBT, resonances and e+e" pairs. 

3. 2-D Gaussian accounting for a mini-jet peak at small A?7, Acj) 

4. cos(A0) term for the dipole on the away-side 

In A-l-A events, this fit function must be augmented with a cos(2A(/)) term to account for an additional component 
that is usually associated with elliptic flow. This simple functional form gives a reasonable description of the main 
features of the data. In Fig. |42l it can be seen that in the more central bins, the long range A77 correlation emerges. 



The parameters of the fit to the mini-jet 2-D Gaussian are shown in Fig. 43 as a function of zy=2Nbin/^part- An 



increase in both the peak amplitude and 77 width, but not the (f> width, above binary scaling (indicated by the dashed 



curves in Fig. 43), occurs for both collision systems and energies at large values of v. It is not yet clear what causes 
the ridge correlation or why there is a rapid onset. This is a particle number correlation, and there is no such rapid 
transition occurring in the equivalent px correlation studies. It would also be of great interest to attempt these studies 



using FID to see if these phenomenon occurs for all baryons and mesons. Table IV shows our estimate of the number 
of events required to perform these studies using non-identified particles. 



TABLE IV: Estimate of events needed to measure the un-triggered ridge correlations It has been estimated that too many 
events are needed to perform this measurement at ^s„„ — 5GeV. 



NN 


11.5 


17.3 


27 


39 


Number of Events 


13 M 


8 M 


6 M 


6 M 



One possible explanation is that this is caused by Glasma flux tube radiation in combination with radial and elliptic 
flow and a flowing sQGP fl50]. A plausible interpretation of the onset of the un-triggered ridge relates the onset to 
collisions passing through the transition scale Qs/ctgiQs) .I'^K- The saturation momentum of partons in the nuclear 
wavefunction, Q^, grows rapidly with nuclear size and energy. Hence as{Qs), the strong interaction strength, is 
correspondingly weak for high energy hadron-hadron collisions and nuclear collisions. 



This could be tested via the energy scan as different ^/s~^ would produce the ridge at different centralities but the 
same transition scale. Fig. 44 and Fig. 45 represent the two cases when the transition scale scales with either energy 
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FIG. 43: Preliminary fit parameters of tlie 2-D mini-jet Gaussian to the pair density correlations for Au+Au and Cu+Cu 
collisions at ^/s~^ = 62 and 200 GeV as a function of i'. From [141| . The vertical lines in the left figure indicate the value of 
f for collisions with impact parameters of b=0. The curves/bands labeled "GLS" indicate the values of the parameters if the 
data scaled with the number of binary collisions. 



density or particle density. The colored curves represent differing collision energies. Below the transition scale, no 
ridge is produced and the curves are interpolated smoothly between the two cases. It can be seen that as the ^/s~^ 
of the collision is lowered, the number of participants required to pass above the transition point rises. In the case 
of energy density scaling no ridge is produced below y/s^^ = 35 GeV, in the particle density case the ridge persists 
down to ^Sjy„ = 13 GeV. These predictions could be tested via the energy scan. 



Ridge Amplitude 




FIG. 44: The un-triggered ridge amplitude as a function 
of Npart when the transition scale is related to the energy 
density. The different colored curves represent different 




Ridge Amplitude 




FIG. 45: The un-triggered ridge amplitude as a function 
of Npart when the transition scale is related to the particle 
density. The different colored curves represent different 



Triggered Correlations 



By triggering on high pT particles and associating them only with the particles in intermediate pt ranges, the jet 
component of the two particle correlations can be selected, while minimizing the soft particle background. In these 
analyses, clear near-side jet peaks are observed |152j as well as the long range A77 correlation 153 . As expected from 
jet fragmentation, the near-side jet correlation yield rises as a function of trigger pt, Fig. |46j but stays approximately 
constant as a function of centrality |154j . The ridge yield on the other hand increases with Npart, Fig. 47 Interestingly, 
while both the jet and ridge yields rise with ^s„„, the ratio Ridge/ Jet is constant, for ^s„„ = 62 and 200 GeV, 
Fig. 48 |154) . It needs to be tested if this is true at other collision energies or merely a coincidence. 
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FIG. 46: The jet correlation yield per trigger on the near- 
side as a function of trigger particle's pr for ^/s~^ = 62 
and 200 GeV Au+Au, Cu+Cu and d-Au collisions. The 
colored curves are predictions from PYTHIA. 



FIG. 47: The number of particles in the near side ridge 
correlation as a function of Npart for ^/s~^ — 200 and 62 
GeV Cu+Cu and Au+Au collisions. 
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FIG. 48: The Ridge/Jet correlation yield ratio as a func- 
tion of Npart for Cu+Cu and Au+Au collisions at .ys„„ 
= 62 and 200 GeV. From IT541. 



FIG. 49: The triggered correlations from SPS compared 
to PHENIX. The away-side distribution does not change. 
Note the scaling of PHENIX by a factor 0.4. 



Both CERES (155 and NA49 |156j have attempted this analysis at the SPS. While they see evidence for a near-side 
jet correlation their away-side peaks are consistent with that arising from momentum conservation only, Fig. [49] Also, 
they do not have the acceptance or statistics to attempt a two dimensional analysis so they can make no statements 
about the presence of the ridge. We therefore propose to pursue this analysis with STAR to further our understanding 
of this phenomenon. 

To estimate the feasibility of these studies we have performed PYTHIA simulations. While PYTHIA does not produce 
the ridge correlation, it is a good representation of fragmentation so we focus on the near-side jet correlation. Since 
the ridge correlation is similar in magnitude to the jet correlation at ^/s^ — 62 and 200 GeV, we determine the event 
rate estimates by assuming that if there are sufficient statistics to see the jet, the ridge should also be measurable, 
i.e. we assume that the ridge/jet ratio remains a constant. PYTHIA simulations were run at various energies and 
the correlations produced via the same techniques used on the real data. Two trigger and associated thresholds were 
used following the cuts of STAR and the lower values of NA49. The cuts used were 1) 3 <pt < 6 GeV/c and 
1.5 <pT < PT GeV/c and 2) 2 <pT < 4 GeV/c and 1 <pT ^"""'^ < pt '"^ GeV/c. The NA49 studies 
are performed with lower thresholds to improve their statistics. This is not possible at the current higher energies of 
STAR because of the increased reach in pT of the soft physics background. The statistics are estimated by requiring 
a better than 50% relative error on the extracted jet correlation yield. The required number of events is shown in 
Table. |Vj This is likely a slight underestimate as the error due to the background subtraction is not included. It can 
be seen that this study is possible down to ^s„„ = 27 GeV when the NA49 approach is used. The ridge yield is 
strongly affected by the V2 subtraction, however as the y/s^ of the collision reduces, so does the measured elliptic 
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flow, so this systematic error is likely to reduce faster than the signal. 

TABLE V: Estimate of events needed to measure high pr jet correlations as a function of y's" 



y^SNN (GeV) 


Number of Events 
2 < pp^ < 4 GeV 
1.0< p^"'"" < pp" 


Number of Events 
3 < p^"" < 6 GeV 
1.5< p^"^"" Kpp" 


200 


0.8 M 


1.6 M 


62 


2 M 


8 M 


39 


4.5 M 


24 M 


27 


8.8 M 


53 M 


20 


23 M 


240 M 



D. Particle Production Studies 



1. Identified Particle Yields, Spectra and Statistical Model Fits 



Since phase diagrams describe only thermalized media, if we are to place a point on the phase diagram we need to 
ascertain that the medium created in heavy-ion collisions is, at least temporarily, in thermal equilibrium. This goal 
continues to be a challenge for the whole heavy-ion community. To this end the bulk/soft physics of integrated yields, 
transverse momentum spectra and particle ratios need to be studied. 

Statistical models have had great success in describing the mid-rapidity, and integrated, particle yields as a function 
of y^s„„ . While successful fits do not prove the validity of the model's interpretation of the parameters as the thermal 
parameters of chemical freeze-out temperature, T^;;, and baryon chemical potential, /i^, should the fits not work 
the system cannot be in thermal equilibrium at hadronization. It is therefore useful to perform these fits before 
embarking on other measures that require the assumption of thermal equilibration to be interpreted. With STAR's 
large acceptance and good PID these measures are achieved to high statistical precision with only a very small event 
count. Previous fits to the data, for example Fig. [78j have shown that from SPS energies upwards the extracted 
Tch is very close to the cross-over temperature, Tc, for a phase transition predicted from lattice QCD of ~ 170-195 
MeV [20 1 1157j . This is one of the indications that the QGP is formed in these collisions. 



2. Strange Particle Yields 



For more than 25 years, it was expected that the production of strange particles in a QGP phase would be enhanced 
with respect to a hadron gas |158l I159j . Also, during hadronisation, the s and bars quarks from the plasma may 
coalesce to form (j) mesons. Production by this process is not suppressed as per the OZI (Okubo-Zweig-Izuka) rule fl60J- 
I163j . This, coupled with large abundances of strange quarks produced in the medium, would again lead to a dramatic 
increase in the production of mesons relative to non-QGP production |158| 1159] . The NA35 and NA49 experiments 
reported jl64l 1165] that in central S-l-S (Pb-|-Pb) colhsions at y^s^ = 19.3 (17.3) GeV the kaon to pion ratio is 
approximately two times higher than in N+N interactions at the same energy per nucleon. However, it was later 
discovered that an even greater enhancement occurs at lower energies, where no QGP is expected |1131I115[I166H174] . 
Interpretation of the experimental data is therefore difficult and additional high precision measurements at lower 
collision energies would be illuminating. 



Another crucial point is the centrality dependence of strangeness production. It was observed by the NA57 and NA49 
collaborations |175H179] that an enhancement is also visible for other strange particle species and that it depends on 



the strange quark content. Fig. 50 shows the comparison of the centrality dependence of S~, A and A yields (|y| < 
0.5) between NA57 and NA49 at ^s„„ = 17.3 (top panel) and ^s„„ =8.8 GeV (bottom panel). For all particle 
species and energies a discrepancy is observed for the more central collisions except for A where a discrepancy is 
also visible in peripheral collisions. Furthermore, the shape of the A centrality dependence at ^s„„ = 17.3 GeV is 
different between NA57 and NA49. 
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FIG. 50: Comparison of the centrality dependence of H , A and A yields (|y| < 0.5) at .ys„„ = 17.3 (top panel) and .^s^ 
=8.8 GeV (bottom panel) between NA57 and NA49. 



Fig. [51] shows the collision energy dependence of mid-rapidity strange particle production from SPS to RHIC measured 
by the NA57 jTTSHTfTj . NA49 [nZlIIED], CERES tlMj and STAR ^[HOIIIHIIIHS] collaborations. A clear discrepancy 
between NA57 and NA49 is visible for all particles. For if ° CERES and NA49 are in good agreement whereas NA57 
measure a slightly higher yield. 

In order to describe the shape and the enhancement in the centrality dependence it is important to consider dif- 
ferent types of models. Statistical models are quite successful at fitting particle yields in central Pb-f Pb/Au-|-Au 
collisions [184] I185j . The centrality dependence of strangeness production is described in terms of a transformation 
from canonical to a grand canonical ensemble in the statistical model |186j . This model uses a correlation volume in 
order to make a comparison to experimental data. A different approach to determine the correlation volume is used 
in the percolation model |187| . In this model, the fireball volume is segmented into smaller sub- volumes (clusters). A 
similar approach is used by the so-called core-corona model jl88j . where an A-l-A collision is assumed to consist of a 
central core in full chemical equilibrium (7<,=1) surrounded by a corona produced by p-|-p type of interactions where 
strangeness is under-saturated. In order to use a hadronic baseline, microscopic models like UrQMD [5S] can be used 
to distinguish between a hadronic and partonic world. So far, discrimination between strange hadron enhancement 
due to dense partonic medium formed in heavy-ion collisions or canonical supression in p-l-p collisions has been highly 
ambiguous, when using the available experimental data on strange hadrons with non-zero net-strangeness. However, 
the measured enhancement of the </> (ss) production (zero net strangeness) in heavy-ion collisions relative to p-t-p col- 



lisions at RHIC, Fig. 52 seems to clearly indicate the formation of a dense partonic medium 189, 190 . Furthermore 
(j) mesons do not follow the strange quark ordering as expected in the canonical picture for the production of other 
strange hadrons |189[ I190j . This seems to rule out canonical suppression effects being the predominant cause of the 
observed enhancement in other strange hadrons in high energy heavy-ion collisions. The BES will allow the further 
testing of the centrality and energy dependence of strangeness enhancement predicted from the various models. 
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FIG. 51: The energy dependence of strange particle and anti-strange particle yields at midrapidity (|y| < 0.5) from SPS to 
RHIC energies. 



3. Mechanisms of particle production at forward rapidity 



At top RHIC energies, two interesting features have been observed at forward rapidities: (a) While at forward 
rapidities particle production scales purely on the number of participants [79| , at mid-rapidity the production follows 
a dependence on both the number of participants and the number of binary collisions [191j . This suggests that there 
are different particle production mechanisms at work as a function of rapidity. The beam energy dependence of 
these scalings needs to be studied to fully understand how particles are produced in heavy-ion collisions, (b) The 
longitudinal scaling of charged particles and n at near-beam rapidities observed at RHIC, left plot in Fig. [53j has 
been found to be beam energy independent, but centrality dependent [1921. The reason for this centrality dependence 
may be attributed to fact that the measured charged particle yields contain significant contributions from baryons 
which vary strongly with the number of participants in the collision. This changing baryon transport plus direct 
contribution from beam protons could be reason for violation of the scaling (as shown in Fig. 53 1. Photon production 
being mostly from tt" decays provides a clean test of the energy and centrality dependence of longitudinal scaling. 



Not many multiplicity detectors will be available at forward rapidity at RHIC during the proposed beam energy 
scan program; data from the PMD will provide crucial information regarding the rapidity dependence of particle 
production. 



4- Light Nuclei Production 



The very small binding energies of nuclei of mass number A imply that these composite objects are formed via 
the coalescence of A nucleons, that are close together in phase space at the time of chemical freeze-out. Their 
production rates and spectra depend sensitively on the interplay of several dynamical aspects of the collisions such 
as nucleon density profiles, temperatures, and collective flow. Volume expansion due to secondary interactions tends 
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FIG. 52: The measured enhancement of various strange particles in A+A compared to p+p at the same beam energy as a 
function of the number of participants. Top paneh For K, (j>, A and H in Au+Au and Cu+Cu colhsions at .ys„„ =200GeV. 
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FIG. 53: Longitudinal scaling of n and protons rapidity distributions for various collision beam energies '193'. Rapidity 
distribution normalized to number of participating nucleon pair is plotted as a function of rapidity shifted by beam rapidity. 



to diminish cluster yields as the overall particle production increases with the total energy ^^s^^ of the collision. 
Counterbalancing effects arise from collective flow which tends to focus nucleons and nuclei in phase space and hence 
increases the production rates of clusters and decreases the apparent source sizes with respect to the total volume 
of the system. Collective motion, temperatures, and densities are also related to entropy production and pressures, 
which are important quantities in the study of phase diagrams. Thus, the study of the rates and spectra of light 
nuclei remains a useful tool for understanding the expansion and correlations of the source. 



37 



A wide variety of physics can be extracted from the study of hght fragment production. 

• A specific ratio of the invariant cross-sections for hght nuclei and nucleons, called "Ba" , can be converted into 
effective source volumes (Va~1/Ba) in a number of model-dependent ways. These effective volumes can be 
compared directly to those obtained from intensity interferometry, e.g. the product of the '"side" and "long" 
radii (i?s2i?i)^/^. The discrepancies in these comparisons are often interesting as they can result from the specific 
model assumptions used to convert the coalescence Ba values into volumes in the presence of strong collective 
flow and the temporal development of the source size (higher mass nuclei freeze-out earlier than lighter nuclei). 

• Light nucleus production rates can be used to estimate particle phase-space densities at freeze-out |194j . Such 
studies have indicated complex relationships between the freeze-out surfaces for pions and nucleons that depends 
on -y/s^ via the overall production rates of pions and nucleons. At the lowest (highest) energies studied at the 
SPS, freezeout is governed by the nucleons (pions) f l94j . 

• The ratio of invariant cross-sections for particles and antiparticles of mass, A, can be used to deduce values of 
the chemical potentials to the temperature (and the square of the "Fugacity"), for example: d/d =exp(-2/i;^/T). 
As the proton chemical potential can be independently measured, the measurement of d/d (and the ratios for 
higher mass nuclei) can provide information on the neutron chemical potential and the n/n ratios, which are 
not generally directly measurable [ 195, 196 . 

• At the SPS (1971 1198] . the initial n/p ratio of 1.54:1 for ^°^Pb evolves towards the equilibration value of ~1:1 in 
the fireball due to strong resonance production. Thus, the measurement of the triton/'^He ratio, for example, can 
be converted using a simple additive procedure to the freeze-out n/p ratio in the source, and hence determine 
the degree of chemical equilibration achieved there. 

• The deuteron to proton ratios can be used to determine the entropy per baryon |199|,I200] via S/A = 3.95-ln(d/p) 
+ (N,/N„). 



Recent studies from the lattice indicate strong effects in the nct-baryon kurtosis (measured via net-protons in 



an experiment) near the critical point (see Section II A 1 1. Similar effects may be visible in the study of the 



fluctuations in the deuteron production rates in the same events. 



The EOS, E877, and E878 experiments at the AGS, the NA44, NA49, and NA52 experiments at the CERN-SPS, and 
the RHIC experiments, have provided a wealth of information on the y/s~^ dependence of light nuclei production 
in A+A collisions. However, a systematic study of the light nucleus production rates and spectra in a RHIC energy 
^s„„ scan may clear up some unresolved aspects of the data from the SPS experiments. These aspects are as follows. 

• The SPS experiments do not in general measure the same kinematic regions near mid-rapidity. The comparison 
of the results thus requires an approximate scaling from one experiment's acceptance to another's, and the 
production rates so compared are not always in perfect agreement |2011 1202| . The Ba ratios measured by NA44 
and NA52 at the top SPS energy differ by a factor of two (a ~3-cr discrepancy) while the NA49 result is in 
between. Only NA49 measured light fragments at SPS beam energies lower than y's^^ = 17.3 GeV. 

• In the NA49 data, the rapidity distributions of '^He nuclei are "concave" versus the rapidity {i.e. indicate a 
shallow minimum at ycM = 0), while the nucleon rapidity spectra are essentially flat in the same rapidity 
region. The observed increase in the '^He formation rates for increasing |ycM| is not yet understood |197l 1198) . 
The observed dependence of the light nucleus production rates on the rapidity in NA44 '203 is a feature that 
disagrees with the model of Ref. 204 , which is otherwise successful in describing direct comparisons of the 
source dimensions inferred from HBT intensity interferometry to those obtained from light nucleus coalescence 
prescriptions. 

• The distinction [205| between the coalescence of A nucleons versus the coalescence of nucleons and lighter 
fragments for the formation of a cluster of mass A is not completely clear. Some analyses imply that deuteron 
formation as an intermediate step to triton formation is required to match the data |206j |207j while others 
disagree |203j . 
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5. (Anti-)Hypernuclei 

A hyper- nucleus is a nucleus which contains at least one hyperon in addition to nucleons. The first hyper- nucleus was 
discovered by Marian Danysz and Jerzy Pniewski in 1952 in a cosmic ray experiment [208 . The smallest and simplest 
hyper- nucleus is hypertriton (3AH), consisting of a A, a proton and a neutron. Hyperons inside a hypernucleus 
contain strangeness quantum number, and hence, provide one more degree of freedom for nuclear spectroscopy than 
the normal nucleus, which consists of protons and neutrons. These hypernuclei also provide an ideal laboratory for 
studying the force between hyperon and nucleon (Y-N interaction), which was otherwise not possible with normal 
nuclei or the traditional hadron-hadron or electron-positron beams. This information is much needed to understand 
the configuration of a neutron star, which depending on the strength of the interaction can alternatively be an object 
of strange quark matter, hyperon star, or a kaon condensate at the core |209j . 

The conventional methods of providing hypernuclei were by cosmic ray interaction, kaon capture or strangeness 
exchange 210 . There are currently several hypernucleus experiments in the major nuclear facilities: MAMIC at 
Mainz and JLab (photo-production), FINUDA at DA(/)NE (stopped kaon beam from e+e" collider), J- PARC (stopped 
kaon from hadron beam), PANDA at FAIR (stopped anti-proton annihilation), and HypHI at FAIR and SPHERE at 
JINR (heavy ion beam). All these experiments require a target of baryon-rich dense nuclear matter to provide the 
nucleons necessary for hypernucleus production, and are incapable of producing anti-hypernuclei: the anti-matter of 
hypernuclei. 

RHIC collisions produce an abundance of particles with a high phase-space density. This environment is thus uniquely 
suited for the production of exotic nuclei and anti-nuclei and anti-hypernuclei via coalescence at the late stages of the 
evolution of the medium produced. The abundances of nuclei and anti-nuclei are similar at the core of the reaction 
zone since the anti-matter and matter are more or less equal at top RHIC energies. This offers the first opportunity for 
discovery of anti-hypernuclei as well as anti- nuclei (A > 3). Since the coalescence process for formation of hypernuclei 
(anti-hypernuclei) requires that nucleons (anti-nucleons) and hyperons (anti-hyperons) are produced in proximity, 
the hypernucleus (anti-hypernucleus) production are sensitive to correlations of the coordinate and momentum phase 
space distributions of nucleons and hyperons. Similarity or equilibrium between these two species is one of the 
signatures of the QGP formation, which would also result in higher hypernucleus (anti-hypernucleus) yields. The 
hypertriton yields can be compared to the yields of '^He and triton which have the same atomic mass number. 
We have found clear evidence for the first-ever observation of a anti-hypernucleus, the anti-hypertriton, as reported 
in |211j . The anti-hypertriton and the hypertriton were reconstructed in the TPC and identified via the secondary 
vertex of hypertriton to (^He -I- tt). We have also been able to measure the anti-hypertriton lifetime, though with a 
large statistical uncertainty (20 to 30%). 

We plan to increase statistics for the signals currently under study, which for example, will result in a more precise 
lifetime measurement for the anti-hypertriton. We plan to carry out a few key measurements of hypernucleus pro- 
duction which will provide crucial information about the hyperon - nucleon correlation. The BES will provide data 
points between AGS and top RHIC energies to establish the trend of the hypertriton/'^He ratio. This can then be 
compared to calculations at a quantitative level. For ^s„„ = 17 (5) GeV, with a penalty factor of 368 (48) |212| and 
a ^He yield of 2x10"'' ( 0.01) and a hypertriton/^He ratio of 0.3 (0.05) the estimated number of min-bias Au-|-Au 
events needed to extract a hypertriton signal at the 5a level is 10-100 M (1-10 M). The above estimates assume that 
the reconstruction efficiency is the same at all energies. In reality, the background and hyperon yields decrease with 
decreasing beam energy while reconstruction efficiency increases due to lower TPC occupancy. 

E. Local Parity Violation in the Strong Interaction 

The observation of local parity violation in the strong interaction in heavy-ion collisions would be an extremely 
important result. The concept that non-central collisions may result in such a violation was first postulated over a 
decade ago |213( I214j . More recently it was suggested that it can result in charge separation relative to the reaction 
plane |215| . Non-central events should produce media with large orbital angular momenta. Since the system has a 
net charge this should result in strong magnetic fields. If the system is also de-confined this may result in P violating 
domains and different number of left and right handed quarks. In the strong magnetic field it will result in preferential 
emission of like sign charged particles from these domains along the angular momentum vector, i.e. in the direction 
normal to the reaction plane. If one could measure the number of positively or negatively charged particles with 
respect to this plane, the P violating domains should result in 
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oc 1 + 2a± sin(0 - "^/rp) + 2vi cos(0 - "iiRp) + 2v2 cos{2{(f) - "^rp)) + ... 



where a± is the asymmetry and thus a measure of the scale of the parity violation. The problem is that while the 
emission direction is preferred to be along the angular momentum vector, the sign of this emission vector is random 
for each domain. This means that {a±) — and since dN±/d(f> can only be measured statistically by experiments 
the proposed signal is lost. What can be measured is the mean pair asymmetry, {aaap), where a and /? indicate all 
possible particle charge sign combinations. 



+ 4)13- ^^Rp)) W [vi.aVi^p - aattp) 



(a+fl-) should show a positive correlation for like-sign pairs and a negative one for unlike-sign pairs [2161 
signals. Fig 
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Such 
200 GeV 



have recently been reported by STAR in non-central Au-|-Au and Cu-|-Cu events at 
where the magnitude of the effect grows with decreasing centrality as predicted |217) . Such a parity violating signal 
only results from the combination of a strong B-field (peripheral events) and de-confinement. This means there is 
likely to be a strong threshold effect that can be sought via the energy scan. Due to the predicted increase in the 
signal as the collision energy decreases (unless the sQGP is no longer formed) it has been estimated that 4 M events 
are sufficient to resolve the signal for Js~l^ = 5-39 GeV. 



Care must be taken however before parity violation is declared in these events because (aaap) is P-even not P-odd and 
so may contain significant contributions from other sources such as resonance decays. Several of these contaminating 
sources have been investigated and none have been able to consistently create signals of the right magnitudes and 
centrality dependence \21S\. The shaded areas in Fig. 54 indicate the magnitudes of the background contributions, 

One way to lend credence to this P-violation measure 
the contributions from the other 



these effects should be smaller at lower collision energies, 
is to observe the same result at different collision energies. By varying 



contaminating sources will occur with different weights. Thus the observation of this potentially P violating signature 
at other energies will help rule out other more traditional explanations of the results. 
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FIG. 54: (cos(<^ci + <t>0 — as a function of centrality for Au-|-Au and Cu-fCu events at ySj^„ = 200 GeV. The shaded 

regions indicate the systematic uncertainties. The thick and dashed curves are HIJING 3-particle correlation calculations for 
Au-|-Au and Cu-|-Cu collisions respectively. Figure from |217| . The uncertainties will be smaller at lower collision energies. 
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III. ANALYSIS OF CURRENT TEST RUNS - AU+AU AT 



9.2 AND 19.6 GEV 



A number of short low energy test runs have been performed at RHIC, the data recorded are used befow to examine 
star's readiness for the BES. In particular data have been taken at ^/s^ = 19.6 and 9.2 GeV in 2001 and 2008 
respectively. The 2008 9.2 GeV run was below injection energy. STAR recorded ~175 k Au+Au collisions at the 
injection energy of y's„„ = 19.6 GeV. From these data 5k events belongs to the top 10% collision centrality and had 
a well reconstructed primary vertex. During the two days of ^/s~^ = 9.2 GeV tests 200 k triggers were recorded 
resulting in ^3.5 k good Au+Au events. A good event is defined as having a reconstructed primary vertex that is 
well within the TPGs acceptance {\z\ <30 cm) and consistent with coming from along the beamline. The results have 
been published in 



Even with these few events the uncorrected minimum bias charge particle tr ansv erse momentum spectrum at 9.2 GeV 



10" 



1(f 



extends to > 3 GeV/c, Fig. 55 and there are indications of an A peak, Fig. 56 
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FIG. 56: The invariant mass distribution of the A for - 
= 9.2 GeV minimumbias Au+Au data. 



While neither of these runs were long enough to improve on the physics measurements reported by the SPS experiments 
many checks have been performed. Some of the results of the analyses of these data are shown briefly below, or have 
been included in the figures and discussions in the previous sections. The resulting px spectra are shown for tt, K 
and (anti)protons in Fig. 57 for the 19.6 GeV and for pious and protons at 9.2 GeV in Fig. 58 |105lll07ll219| . Similar 
results exist for the charged kaons. We measure 82% of the produced tt at mid-rapidity and 75% of the protons, 
the kaon reach is not as complete, covering only 47% of the yield, however this will be greatly improved with the 
inclusion of the ToF. HBT radii have been extracted from the y/s~^ = 9.2 GeV Au+Au collisions for tt in the 0-30% 
centrality bin, Fig. 60 The fits yield results that correspond with the ^s„„ systematics from AGS to RHIC shown 
in Fig. 61 |220j . No anomalous jump is seen in the Rout/R-side ratio, although the peak could be narrow and located 
between y/s^ =20 and 60 GeV. The feasibility of such a measurement is indicated in Fig. 59 where statistically 



relevant identified V2 measurements have been made with only ^3 k events. To date all of STARs results fit well into 
the measured systematics reported by the SPS, indicating that STAR is ready for full length data taking. 



Running at sub-injection energies requires a different harmonic to top energy running and is challenging to both the 
experiments and the collider. For instance the particular choice of ^s„„ = 9.2 GeV did not allow the beams to be 
cogged at STAR and PHENIX simultaneously. However, data were recorded and this cogging issue will be resolved 
for the BES by an appropriate selection of collision energies. The data taken during this short run are used below 
to demonstrate that the event rate projections for the BES are correct, and that STAR's triggering, and event and 
vertex reconstruction capabilities, designed for top energies, are applicable at these low y/s~^. 

During this test beam an event rate at STAR of ~ 1 Hz was observed. This is a much lower rate than that expected 
during the actual energy scan due to a number of factors. The two major ones were: a) since this was a test run 
only 56 bunches, out of a maximum of 120, were used and b) the intensity of each bunch was less than the maximum 
possible by a factor 3-6. A further gain in measured luminosity could be obtained by running in continuous injection 
mode. All these factors combined suggest that during the actual energy scan an event rate of >5 Hz can be expected. 



Two y/s^ = 9.2 GeV Au+Au events are shown in side and end view in Fig. 62 and Fig. 63 The collision vertex 



is clearly defined and the mean 77 of the charged tracks is approximately zero in both cases, strongly indicating that 
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FIG. 57: The identified particle spectra as a function of centrality for Au+Au collisions at ^s„„ =19.6 GeV. From 



these are beam-beam collisions and not beam-gas or beam-beampipe. 
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FIG. 58: The corrected a) tt and b) p pr spectrum for 
y'Sj^j^ = 9.2 GeV Au+Au data as a function of centrality 
at midrapidity {\y\ <0.5). From 



FIG. 59: The charged particle V2 at 
function of pT- From [TII51[^ . 



=9.2 GeV as a 



Many of the events recorded were beam-beampipe collisions due to the large diameter of the 9.2 GeV beam. These 
events, Fig. [64j are easily rejected off-line via cuts on the radial position of the primary vertex. In the left plot the 
beam spot due to beam-beam collisions is clearly seen at x ~ 0.5 cm and y ^ -0.5 cm. The beampipe, radius 5 
cm, is also illuminated due to beam-beampipe events. The right hand plot of Fig. |64] shows the location in z of the 
beam-beampipe collisions. STAR's beampipe is made of beryllium in the \z\ < 50 cm region and changes to aluminum 
in the \z\ > 50 cm, outside of the nominal collision diamond. The increase in the number of events occurring in the 
large z range is due to the higher cross-section of Au on aluminum. Since the event rate will be low it has been 
decided that no attempt will be made to reject these online. 



IV. THE P+P BASELINE 



For certain measurements in the A-l-A energy scan (^s„„ =5 — 60 GeV) it is important to know well the baseline NN 
cross section and invariant spectra to use as references. Previous relevant measurements come from many different 
accelerators, AGS, PS, SPS, FNAL and ISR spanning 40 years of experiment. The quality of these measurements are 
in general poor, particularly at higher px due to limited statistics, but this may be all that we will have access to, since 
performing these measurements at RHIC is not feasible in the first years of BES program. In addition, agreement 
between different experiments is not perfect in several cases with differences up to 50% at high pr- A number of 
articles have discussed the p-l-p systematics from available data in particular for high-pT particles in regard to the 
RHIC data. Most data are for pions, both charged and neutral. Data are much sparser for other hadrons. 



For fixed target energies relative to y^s~^ we tabulate the relation between the beam momentum and the ^/s~^ for 
a number of energies where experiments have been performed in fixed target experiments, Table VI In the following 
sections we discuss briefly relevant pion data from p+p collisions at lower energies, particularly related to the higher 
Pt ranges. 
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FIG. 60: The HBT projections of the 3-dimensional cor- 
relation functions for negative vr data from 0-60% central 
Au-|-Au collisions at .^s„„ = 9.2 GeV. Fits shown are the 
Bowler-Sinyukov functions [22l[223]. From [T05l[22T] . 



FIG. 61: The HBT systematics as a function of .Js^ 
Low energy data from STAR added to plot from [220] . 
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12 


4.55 


CERN 


24 
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CERN 
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17.3 
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200 


19.32 


CERN,FNAL 


400 


27.36 


CERN,FNAL 


800 


38.8 


FNAL 



TABLE VI: Relation between P;ai, and cm. energy in p-f-p collisions, where pp data exist. 



A. Energies below top SPS 



The very low energies were studies at AGS and CERN in the 1970's. In particular data were obtained at 12 and 24 
GeV incident energy. Production at these low energies has no hard components, but arises essentially from various 
resonance productions resulting in an approximate m^-scaled spectra up to the kinematic limit. This same feature is 
seen in A-l-A collisions, see e.g. Ref. |224j with preliminary data from Si-|-A at 14.6 A GeV. The main p-|-p data comes 
from Ref. 22"i and for higher data from Ref. |226) . In the energy range of 24-158 GeV fixed target there is really 
no available data and, as will be discussed, extrapolating between the measured data may be difficult and imprecise. 



B. SPS, FNAL and ISR energies 

A fairly extensive data-set for pion production from NA49 [227] covers a large xp range. The data at mid-rapidity 
extend to pr of 2.1 GeV/c. A summary of data around ^s„„ = 22 GeV was made by Arleo and d'Enterria j228, in 
order to parametrize the cross section for use in the comparison to RHIC data at y/s^ = 22A GeV (Cu-|-Cu). They 
made small energy dependent corrections based on pQCD to the available data between 21 and 23 GeV and gave a 
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parameterization that will be used below. The parameterization is given by 

f{pT)^p^\l + {pTlpi)Y^*[l-XTr, 

with po = 176.3, pi==2.38, P2==-16.13, Ps. = 6.94, and xVndf = 208.2/190. 

At FNAL high momentum 7r+and 7r~data from p+p collisions were obtained by the E605 experiment at fixed target 
energies of 400 and 800 GeV with significant integrated luminosity, 436 p&~^ and 615 respectively |229j . The 

y/s~^ =38.8 GeV data (800 GeV fixed) probes pt out to 10 GeV/c, i.e. x^ of ~ 0.5. The errors on the data are 
typically fO% statistical and a similar value for the systematic uncertainty. This data-set is fairly consistent with a 
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FIG. 64: Left: The radial position of the reconstructed primary vertex in the Au+Au ys^^^ = 9.2 GeV test run. Right: The 
z position of the primary vertex for those events with a radial primary vertex position indicating they occurred in the beam 
pipe, R > 2 cm. The beampipe is constructed from beryllium for \z\ <50 cm and aluminum outside of this range. 



measurement by the earlier data from the Chicago-Princeton collaboration |230j . albeit there are differences in the 
order of 30%. The data at ^s„„ =27.4 GeV has smaller pT coverage. 

The ISR tt" data has often been summarized emphasizing the a:T-scaling properties, as has the fact that while the 
within experiment energy dependence (like CCOR) are well established, overall 30% disagreements on cross sections 
are quite common (see e.g. Ref. |231j ). This kind of disagreement also exists between the ISR and FNAL experiments, 
even though the general trends are in good agreement. 
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FIG. 65: Invariant cross-sections for tt , vr", and vr for var- 
ious cm. energies as indicated in the legend. The data are 
from reference given in the text. 



FIG. 66: Scaled invariant cross sections for several energies. 
The full drawn curve is the best fit to the 22.4 GeV data. 
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C. RHIC 62.4 GeV 



After the first RHIC A+A run at 62.4 GeV (Run 4) a compilation of p+p data for that energy was made by d'Enterria 
[232j and used for the first i^AA results. The different measurements from the ISR show significant differences from 
the mean cross sections (vs pr) so a dedicated p+p run at 62.4 GeV at RHIC was warranted and carried out in 2006. 
The ISR data agree reasonably well with data from PHENIX [233] and unpublished BRAHMS data on 7r+ and tt^. 
It noted that PHENIX changed the energy scale of the data from the CCOR experiment |234) . 



D. Energy dependence of high-pT cross sections 



In Fig. [65] we show some selected data for pion production from AGS to ISR energies plotted vs. pr- The typical 
decrease of the hard scattering cross section as the energy is lowered is observed. Ref. f226] presents an overall function 
that describes the available data as of 1978. It gives a good overall description of data from AGS to ISR energies. We 
refer to the reference article for the functional form, which is quite complicated. Since its publication additional data 
has been reported and for the specific formulae to be used as a reliable interpolation between energies ^/s^ = 5 — 30 
GeV it will have to be revisited and refit. However, it is expected to be well suited for a first order estimate of cross 
section in this energy regime. 

In Fig. [66] we present the data from Fig. [65] in terms of the variable xt as has been customary. STAR and PHENIX 
demonstrated in |233l 1235] that the 62 and 200 GeV data falls on a common curve at large xt when the cross sections 
are scaled by y/s~^ " with n=6.5 ± 08 and 6.38 respectively. At the lower ISR and SPS energies an increased n 



factor is needed to describe the data well. In Fig. 66 the cross sections have been scaled by \/22.4 with n=8.2 as 
indicated by the CCOR data |234j . The very low energy data do not follow this trend, but correspond to an even 
higher factor of n indicative of non-perturbative effects dominating at lower energies. A value of n^ 11 would allow 
scaling at high values of Xt > 0.4 from the two energies. We conclude that for Xt > 0.15 good scaling with energy is 



observed; thus the scaling function from Ref. '228^ (solid curve in Fig. 661 is a usable description. Of course it cannot 
be used at lower xt- 



In addition the NA61 (SHINE) experiment is approved to run at the CERN SPS 236 . In 2009 they collected between 
2-6 X 10^ p+p events at = 6.4, 7.8, 9.0, 12.7 and 17.8 GeV. These new data should supplement those currently 

available, providing high statistics data at very similar collision energies to those planned for the RHIC BES. 

In conclusion existing p+p data and parameterizations for ^^s^^ > 17 GeV gives an adequate, albeit not perfect, 
description. At lower energies the higher px cross sections are not well known and subject to large interpolation 
errors. This may not be a big concern since at RHIC with the low luminosities at the low energies of the first BES 
this momentum range will not be accessible. 



V. THE STAR DETECTOR 



The STAR collider geometry immediately provides advantages in comparison to the SPS fixed target experiments. 
These advantages are two fold [237] . The first is that the detector occupancy at mid-rapidity increases much faster 
as a function of ^/s~^ in a fixed target experiment than in a collider environment, as shown in Fig. 67 The second 



is that the particle acceptances are dependent on the beam energy in fixed target experiments, e.g. Fig. |68[ whereas 
they are constant at colliders and similar, at mid-rapidity, for all particle types, Fig. [69] and Fig. [70] The reduction in 
track density means that for a given collision energy a detector in a collider environment has less problems with charge 
sharing of hits, and track merging resulting in a better reconstruction efficiency and cleaner particle identification. 
The constant acceptance also means that there is better control of the systematics of the measurements and that a 
number of the uncertainties cancel when comparisons at different ^ys~^ are made. This is especially important for 
the energy scan as the identification of the CP is likely to be made by observing how various measures alter as a 
fimction of y^s^. 
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Njs;^ (GeV) 



FIG. 67: The mid-rapidity track density at 1 m radius in a collider setup compared to that of a fixed target as a function of 
^ys^j^. Figure from |237| . 
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FIG. 68: The acceptance for n and K used in the fluctuation analyses as a function of pr and rapidity for NA49 at SPS fixed 
target collision energies of 20 and 40 GeV. Figure from [237] . The acceptance is greater for the inclusive spectra analyses. 



A. STAR'S Subsystems 

This section briefly describes the various subsystems of STAR relevant to the BES, Fig.[7T] Of particular interest to 
the CP search is the installation of the full barrel ToF which was completed during the 2009 summer shutdown. For 
a full description of STAR see |238j and references therein. 
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FIG. 69: The acceptance for tt as a function of pr and 
rapidity for STAR for = 9.2 GeV. 



FIG. 70: The acceptance for K as a function of pr and 
rapidity for STAR for ys];;^ = 9.2 GeV. 




FIG. 71: Schematic of the STAR detector showing aU the subsystems, including those currently being proposed and undergoing 
research and development studies. 



Trigger Detectors 



We propose to use the Beam-Beam Counters (BBCs) for triggering during the low energy runs. The collision rates 
at these low ^ys'^ are only a few Hz so all events can be recorded and no selective trigger is required. Therefore it 
is proposed to use only a minimum bias trigger, i.e. that of a coincidence in the two BBCs. 



The BBCs, shown schematically in Fig. [72j are positioned at ± 3.5 m from the center of the TPC along the beam 
direction. The small inner tiles, shown in blue in Fig. 72 cover 3.8 < |?7| < 5.2 and can all be inscribed within a 
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circle of 9.64 cm. The outer (red) tiles cover 2.2 < \r]\ < 3.8. The RHIC beam passes through the center of the BBC 
(labeled B in Fig. 72 These annular scintillator detectors are sensitive to charged tracks down to the single minimum 
ionizing particle (MIP). Even at the lowest collision energies the number of produced particles in heavy- ion collisions 
is larger than in p+p collisions at y/s~^ = 200 GeV where they perform well. Tabic VII shows estimations of the 



number of particles hitting the inner and outer sections of the BBCs as a function of the mteraction's centrality for 
Au+Au collisions at ^s„„ = 5 and 8.75 GeV. 



STAR Beam-Beam Counter Schematic 

Front View 




Version 4(16/01 -2 
Updated 2/25/02 



FIG. 72: A schematic diagram of one of the BBCs. 



TABLE VII: Estimation of number of charged particles in the BBCs in Au+Au collisions at .ySj^„ — 5 and 8.75 GeV for 
various impact parameters. 



Impact Parameter 


^s;;;^ = 5 GeV 


^/s^ = 5 GeV 


^s];;^ = 8.75 GeV 


= 8.75 GeV 


(fm) 


BBC Inner 


BBC Outer 


BBC Inner 


BBC Outer 


0<b<3 


5 


27 


12 


54 


3<b<6 


11 


30 


21 


57 


6<b<9 


22 


35 


39 


40 


b>9 


44 


30 


66 


8 



The BBC's will also be used to calculate the first order event plane, although with reduced resolution to that calculated 
from the FTPCs. More details are discussed below. 



2. Time Projection Chamber 

The main tracking detector in STAR remains the Time Projection Chamber (TPC). The active dimensions of the 
TPC radially are from 50 cm to 200 cm and it is 400 cm in length. This gives the TPC full azimuthal coverage 
and a uniform acceptance within |?7| < 1. It is placed within a uniform solenoidal magnetic field which allows for 
the determination of the momenta of charged tracks. When the magnet is set to its maximum value of 0.5 T the 
transverse momentum resolution is Apr/pr ^ 0.012pT (pt in GeV/c). 
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3. Time of Flight 



The Multi-Gap Resistive Plate Chamber Time of Flight (MRPC ToF) [2391 will completely surround the outer 
radius of the TPC, |7y|<0.9, 0<(j)<2Tr. It consists of 23K channels from 120 modules. The first MRPC tray was 
installed in Run 3 and it was full installed for Run 10. The ToF has an intrinsic stop time resolution of 80ps. After 
folding in STAR'S start signal a total resolution of at least 120 ps will be achieved. This 120 ps total resolution is 
for the lowest ^s„„ values and will decrease significantly for higher energies. This gives the ToF an excellent mass 
resolution, see for example Fig. [73} 



4- Forward Time Projection Chamber 



The Forward Time Projection Chambers (FTPCs) are two radial drift TPCs which cover the region 2.5< \r]\ <4.0. 
They have a 2-track resolution of 12 mm and a momentum resolution of 12-15%. Unfortunately due to the small 
number of hit points, there are only 10 recorded per track, in the FTPCs no PID is possible. It is planned to use the 
FTPCs as the main detectors used in the determination of the event plane, see later. 



5. Photon Multiplicity Detector 



The Photon Multiplicity Detector (PMD) |241j measures the inclusive number of photons (dominantly from 7r° 
decays) produced in the pseudorapidity regions of -2.3 to -3.8, with full azimuthal coverage. It also provides the 
spatial distribution of the photons in (77, (f>) phase space. 

The PMD therefore provides a unique opportunity of addressing physics at a higher baryon chemical potential (/is) 
compared to midrapidity measurements at a fixed beam energy. It also allows us to continue the studies of various 
scalings with Npart, y Ybean-n Started by the now decommissioned PHOBOS and BRAHMS experiments. In con- 
junction with the FTPCS (both of which are needed), it will provide a unique opportunity to look for a Disoriented 
Chiral Condensate and hence chiral phase transition. Measurements of vi and V2 at forward rapidity for photons and 
neutral pions can be made. 



B. Particle Identification 



With completion of the full barrel ToF STAR will have excellent mid-rapidity PID capabilities. Good quality track- 
by-track particle identification is necessary for our proposed fluctuation measurements. 



Charged particle identification will be performed via a combination of measured ionization in the TPC (dE/dx) and its 
ToF. The addition of the ToF allows for the identification of particles where their dE/dx as a function of momentum 
merge. Using these two techniques STAR will be able to resolve K and tt up to momenta of 1.6 GeV/c and protons 

shows 
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up to 3 GeV/c. Fig. 73 shows the measured 1//3 of the ToF in d-Au collisions at ^/s~^ = 200 GeV and Fig. 
the measured dE/dx of tracks in the TPC, the tt, K, p and electron bands are clear. The colored curves indicated the 
± 8% deviation from dE/dx calculations using the Bichsel parameterization |242j . Using ionization measurements in 
the relativistic region of the TPC and the data from the ToF also allows the statistical identification of protons and 
TT up to 12 GeV/c 

In addition STAR can use secondary vertexing to identify K^, A, A, S^, S^, il~ , and Q~ over the whole pr range, 
and, via invariant mass mixed event techniques, a variety of resonances such as the </> and K* and p ( e.g. |244j and 
references therein). 



C. Event Plane Resolution 



STAR has excellent event plane resolution even at low y/s 



Fig. 75 shows the measured second order event plane 
= 9.2 GeV. Also shown for comparison 



resolution using the TPC as a function of centrality for Au-t-Au events at y/s~^ 
is the event plane resolution for NA49 at y/s^ =8.75 GeV ^91j. It can be seen that STAR's larger acceptance results 
in a much superior resolution. While the V2 event plane can be calculated using the TPC, for vi measurements an 
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FIG. 73: Particle identification using the ToF in d-Au colli- 
sions at ySj^„ — 200GeV. tt, K, p, and electron are clearly 
identifiable. The inset shows the projection of the ToF 
mass^ measurements in the 1.2 <pt < 1.4 GeV/c range. 



FIG. 74: dE/dx measurements using the TPC. The curves 
represent the ±8% deviation from dE/dx calculations using 
the Bichsel parameterization for each particle. 



independent determination of the event plane is necessary, this can be done by utilizing the FTPCs or with slightly 
less resolution the BBCs. Calculating the event plane in this away also minimizes several auto-correlation effects. 



c 
o 

o 

CD 

o 
O 



0.6 



0.4 



0.2 



2 order eventplane resolution 

STAR (9.2 GeV) 







20 40 

Centrality (%) 



60 



5Sl00r 

§90^ 
& - 
80 r 

70 r 

60 r 

50 r 

40 r 

30 r 

20 r 

lOr 















A 


RQMD:5GeV 














T 


























































, * 




















~«- • ■ 
















A 




































▼ 


r 
















































.1. 






( " 2 


t e 


f 


i 1 


1 


2 14 



b in fm 



FIG. 75: Second order event plane resolutions for Au-|-Au 
events at ^/s~^ =9.2 GeV compared to that of NA49 at 
^s];;^ =8.75 GeV [5T]. 



FIG. 76: Estimated first order event plane resolutions as a 
function of impact parameter from RQMD simulations of 
Au-I-Au events at ,ys^^ = 5 and 10 GeV using the BBCs. 
Similar results are obtained for all .^s„„ . See text for 
details. 



The first order event plane resolution of the BBCs is shown in Fig. [76] as a function of impact parameter for RQMD 
simulations of Au+Au events at y/s^ — 5 and 10 GeV. For this simulation the BBC was fully instrumented. For the 
2008 test run at y^s^ = 9.2 GeV only the inner tiles were present so the resolution was degraded, (for details of the 
BBC please see section V A 1 ). The resolution for this dataset was estimated to be ^ 25% which is in good agreement 
to the resolution obtained from RQMD when the BBCs were simulated without the outer tiles. This confirms that 
the simulations are a good representation of the data. The resolution is similar for all collision energies discussed in 
this proposal. 
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D. Acceptance, Efficiency, and Particle Production Rates. 

In order to make realistic estimates of the number of events needed to perform the analyses described in the previous 
sections it is important to be able to estimate the acceptances and tracking efficiencies of the various particles as well 
as their expected yields. 

The correction factors (acceptance*efficiency) used for tt, K and protons are estimated from those calculated for the 
^/s^ = 19.6 GeV test run. For the (p, K^, A, and A results have not yet been obtained for the y/s„„ = 19.6 GeV 
data. Therefore the cj) correction factors obtained for the most peripheral (70-80%) Au+Au at y/s~^ = 200 GeV 
data were used, and for the A, and A the Cu+Cu at ^/s^ — 200 GeV values were used. For all particles that 
decay the branching ratio to their charged daughters are included in the correction factors. All curves were fit to the 
parameterization 



Eff = a X exp[{-b/pTY 



where a, b and c are free parameters. While there is a correction factor dependence on occupancy it is small over 
these energy ranges and therefore neglected in these approximations. The resulting correction factor curves are shown 
At the lowest 



in Fig 77 



\/Snn proposed the efficiencies may be slightly lower than those assumed here due to the 



expectea slightly worse timing of the ToF, see section V A3 




FIG. 77: The mid-rapidity correction factors as a function of pt assumed for the BES. 



The estimated acceptances and efficiencies can be folded with the expected px spectra to determine the raw recon- 
structed particle yields as a function of ^/s^ and centrality. To estimate the rates and (pr) at mid-rapidity for 
several .^s„„ energies and 'Npart a statistical model can be used together with the following formulae, from Ref. |245j . 



For a given y/s]. 



^ 1.308 

1 + 0.273^^ ^ ' 

T = 0.166 - - {)mZ^i% (2) 
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This gives a very good phenomenological parameterization of the current data as shown in Fig. 78 |245j . Having 
calculated T and hb the statistical model can be used to calculate each particle/Nc^ ratio. Nc^, is then calculated via 
a phenomenological parameterization of Glauber model calculations of the number of binary collisions, N^in, and the 
number of participants, Npart as a function of ^/s~^. 



Nun = (0.314 + 8.7e-4yW- 1.8e-6y^^)7V 



part 



(3) 



N,h = (0.5933?n(ys]^) - 0.4153) x ((1 - 0.11)7Vp„rt/2 + O.llA^fc, 



(4) 



Given the particle/Nch ratio and Nc/i one can then estimate the mid-rapidity particle yields. A Blast- Wave parameter- 
ization is used to calculated the {pt) where Tfci„, the kinetic freeze-out temperature can be assumed to be independent 
of ^/s^ and set to 0.1 GeV/c while the mean transverse velocity can be estimated from: 

(/3) = (0.388 + 0.186 ln(0.7928iVpart + 0M29N^^^t - 3. 316006'^ N^^^^ + A.01681e-^N^^^t)/Npart) 

(0.05727 In(^A^) + 0.2933) /0.592 (5) 
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FIG. 78: The calculated chemical freeze-out temperature, T, and chemical potential, /^s, as functions of -v/s^vw 
are the phenomenological parameterizations used in our estimates of the particles yields. Figures from |245| 



The curves 



VI. SUMMARY AND THE RUN PLAN 



For the first BES we propose to run a wide range of energies, from ^s„„ = 5-39 GeV, to truly scan over the available 
collision energies. We have selected lower y/s~^ values that give the greatest discovery potential for the CP as well 
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as higher energies to cover the current gap between the SPS and RHIC. This wih ahow us to explore the questions 
about "turn on/ofF' effects of signatures of partonic media such as the high and low px ridges and constituent quark 
scaling of the elliptic flow. Table |VIII| summarizes the measurements that drive the number of events we request at 
each collision energy. 



TABLE VIII: Estimate of number of events needed at each collision energy. The most statistics hungry measurements are 
listed. 



y/s^^ / Number of events 


5 GeV 


7.7 GeV 


11.5 GeV 


17.3 GeV 


27 GeV 


39 GeV 


V2 up to ~ 1.5 GeV/c 


0.3 M 


0.2 M 


0.1 M 


0.1 


0.1 M 


0.1 M 


Vl 


0.5 M 


0.5 M 


0.5 M 


0.5 


0.5 M 


0.5 M 


Azimuthally sensitive HBT 


4 M 


4 M 


3.5 M 


3.5 


3 M 


3 M 


PID fluctuations 


1 M 


1 M 


1 M 


1 M 


1 M 


1 M 


Net proton kurtosis 


4 M 


4 M 


4 M 


4 M 


4 M 


4 M 


Differential corr. and flue, vs centrality 


4 M 


4 M 


4 M 


4 M 


4 M 


4 M 


Nq scaling tt, K, p, A to ttlt — mo/Nq ~ 2GeV/c 


N/A 


6 M (marginal) 


5 M 


4.5 M 


4.5 M 


4.5 M 


Nq scaling (j), Q, to pr/N, ~ 2GeV/c 


N/A 


56 M 


25 M 


18 M 


13 M 


12 M 


Rep to pT ~ 4.5 GeV/c (17.3), 5.5 GeV/c (27), 6 GeV/c (39) 


N/A 


N/A 


N/A 


15 M 


33 M 


24 M 


Un-triggered ridge 


N/A 


27 M 


13 M 


8 M 


6 M 


6 M 


Local Parity Violation 


N/A 


4 M 


4 M 


4 M 


4 M 


44 M 



Rate estimates are based on a number of ingredients. Most important are the data from the Au+Au beam test at 
y/s^ = 9.2 GeV taken during Run8. In particular, these data allow a determination of the actual rate of events 
on which STAR can trigger and that lie within a range of vertex positions for which STAR can perform the desired 
physics analyses. Thus, rates can be estimated as the number of good physics events with less uncertainty of the 
usable fraction of the total luminosity. The test run was very short with the emphasis on demonstrating capability 
rather than tuning for the highest possible event rate. The BNL Collider Accelerator Division (CAD) staff have 
indicated a high degree of confidence that these rates can be increased by a factor of about 6 via improvements 
in injection efficiency and increasing the number of bunches in the machine. Additional tuning is likely to provide 
further incremental rate improvements. Another option being explored is continuous injection to fill the bunches and 
extend the beam lifetime. This is expected to increase the integrated rate by a factor of about 2. Very significant rate 
increases may be possible using electron cooling. Until the details of the beam loss mechanisms at these low energies 
are determined, the exact gain that could be achieved is unknown but factors of about an order of magnitude are 
considered very possible. To be very conservative, none of these further possible enhancements are included in the 
estimates presented below. 

Again based on consultation with CAD 246j . it is estimated that the rates will scale up and down from the 9.2 GeV 
values by 7'^ up to the injection energy of and by roughly 7^ above that. To cross-check these estimates, there are 
actual data rates taken by STAR during the early injection-energy run and also measurements of the injection energy 
luminosity under current operating conditions. The values typically quoted are total energy, so the 7 of each beam 
is found by dividing by 2 and then by the average nucleon mass for a gold nucleus. These factors are independent 
of energy and therefore, the rate below injection energy simply scales with the cube of the beam energy. Using 
the measured values from the test run increased by the conservative values listed above, the rates are the following: 
Below injection energy: 6.5x10^"^ x (.^s„„) events/sec At injection energy: 48 events/sec Above injection energy: 

12.6x10^^ X (^Sjvn)^ events/sec Note that, as mentioned above, these are estimates of the actual rate of events for 
which STAR can perform all of the desired physics analyses and are therefore not estimates of the total luminosity 
available from the accelerator. These estimates imply the following times required for 1 million events assuming 10 
hours per day of beam. 5 GeV 0.8 evts/sec 35 days 10 GeV 6.4 evts/sec 4.3 days 20 GeV 50 evts/sec 0.6 days 40 GeV 
200 evts/sec 0.14 days 



A suggested run plan is shown in Table IX The length of the run at each energy is determined by the minimum required 
number of events to perform the detailed measurements discussed in the previous sections. The lower beam energies 
are specifically chosen to map out region around the "horn" in the K/tt ratio observed by the SPS experiments |247) . 
All the selected energies allow collisions at both STAR and PHENIX; the lower y/s~^ values are also mapped as 
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closely as possible to those already provided at the SPS. 



TABLE IX: A run plan for STAR assuming a 10 hour day. 



Vsjvjv (CeVJ 


MS (MeV) 


Rate (Hz) 


Event Count 


Run Time (days) 


5.0 


550 


0.8 


Beam development 


7 


7.7 


410 


3 


4M 


36 


11.5 


300 


10 


4M 


15 


17.3 


229 


33 


15M 


13 


27 


151 


92 


33M 


10 


39 


112 


190 


24M 


4 



After analysis of the first run period is mature we propose that a second scan is performed focussed more specifically 
on a few collisions energies. These energies and physics topics will be chosen to explore in more depth the most 
interesting regions found via this first scan. We propose two runs since the second run, more data at fewer energies, 
can take advantage of further luminosity upgrades proposed by the RHIC accelerator division. 

In summary there is great discovery potential for a Beam Energy Scan at RHIC. The most exciting result would be the 
identification of a Critical Point. Equally interesting would be the proof that at lower y/s^ energies the transition 
to a QGP occurs via a first order phase transition. Either of these results would enable us to make precise entries on 
the Phase Diagram of QCD matter. Currently the most promising predictions for identifying these phenomena are 
critical point fluctuations of cither pT or conserved quantum numbers (baryon or strangeness), pair correlations or 
non-monatonic deviations in the energy dependence of flow characteristics. A scan will also provide results that would 
enable a more precise determination of the region where sQGP medium effects, such as constituent quark scaling of 
elliptic flow, (dis) appear. 

Finally we have demonstrated that not only is STAR committed to performing a Beam Energy Scan but that the 
detector system is ideally suited for such an explorative endeavor. With the run plans detailed above we can make 

significant improvements not only to the existing SPS data but also plan to extend the investigations into regimes 
not yet probed by any experiments, and/or to include unexplored physics observations. 
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